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I . INTRODUCTION 


It is the objective of this work to describe the dynamic 
behavior of footings and to develop practical methods for 
dynamic analysis. The Investigations are limited to the 
simple case of a rigid circular footing resting on a plane 
soil surface and excited by a vertical time-dependent force 
in the axis of symmetry. 

In order to arrive at a model which can be analysed 
mathematically it is assumed that the subsoil can be consid- 
ered as a perfectly elastic, isotropic and homogeneous half- 
space and that only normal stresses are transferred at the 
interface between footing and soil. 

Half-space models have been used by several investiga- 
tors, notably Relssner (l) who, in 1936, developed a method 
for finding the steady- state surface displacements of an 
elastic half space for any given surface pressure distribution 
with axial symmetry. Through the work of other investigators, 
such as Sung (2), Bycroft (3) and Hsleh (4), the art has 
since been developed to the point where it is possible to 
calculate the steady state response of a rigid footing-soil 
system in the low frequency range. Fortunately, this range 
includes the operating frequencies of most machinery and the 
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above solutions are therefore useful for many practical r-Jr- 
poses. This was clearly shown In a paper by Rlchart ( 5 )., 
who also considered the rocking and sliding modes of vl ora- 
tion . 

In the present work it is shown how the frequency range 
can be extended to Include all frequencies of steady-state 
motion and also how this complete steady state solution can 
be used to calculate the transient response caused by an 
arbitrary pulse loading. The analytical solutions compare 
well with test results by Pry (6), Chae (j), and Drnevich., 

Hall and Richart (8). 

The above system differs from most classical dynamic 
systems in that it is of infinite dimensions. This geomet- 
rical peculiarity causes an apparent loss of energy by wave 
propagation into the half space. It will be shown in this 
work that the effect of the energy loss is comparable to that 
of the dashpot in a simple damped oscillator. This observa- 
tion has led the author to the adoption of a simplified 
sprlng-dashpot analog for practical calculations. The new 
analog differs from previously suggested analogs in that its 
components (mass, dashpot, spring) are Independent of the 
frequency of the exciting force. The analog can therefore 
be used to predict the response to dynamic pulses of arbi- 
trary shape by means of the classical methods for a damped 
oscillator (phase-plane methods, numerical Integrations, etc.). 



II . STEADY- STATE MOTION 


General Theory 

Before focusing the attention on the elastic half-space 
model it Is appropriate to study a larger class of dynamic 
systems . 

Consider a linear dynamic system S which is excited by 

a periodic vertical force P(t) of frequency u) and amplitude 

P . That is: 
o 

^ ( 1 ) 

Furthermore j let it be assumed 
that P attacks at a point 0 
which is such that the dis- 
placement S of point 0 is 
vertical at all times. The 
class of dynamic systems thus 
defined clearly Includes all 
systems which are^ axially symmetric about a vertical axis 
through 0, and the system.s which will be considered in the 
following are of this type. 

The system S may or may not exhibit viscous damping and 



Pig. 1 
System S 
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it may be of finite or infinite dimensions. 

In the case of steady- state motion it is knovm that all 
forces and displacements are harmonic with the circular fre- 
quency c<J> and proportional to the amplitude of the excit- 
ing force. The displacement 6 can therefore be xvritten in 
the form: 






( 2 ) 


where k is a quantity having the dimension force/length and 
P is a dimensionless function. 

The quantity k will be called the spring constant and 
can usually be put equal to the static spring constant of 
the system S. There -are however members of the above clar.s 
of systems for which this is not possible (l.e., if S con- 
sists of Just a dashpot) and in this case a different choice 
must be made. 

The time-independent, complex function F = P,, + iP^ 
is in general a function of the frequency 60 and the prop- 
erties of the system S. It will, in the following, l )0 re- 
ferred to as the displacement function *for the system S. 

It should be noted that if k is the static spring con- 
stant then we must have P = p^ = 1 for cu =0. 


The displacement function is proportional to the quantity 
compliance used by authors favoring electro-mechanical 
analogies . 
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It is of Interest to Investigate if there exists a 
simple damped oscillator which as far as the displacement 

S is concerned, can be used 
as an analog for the system S, 

In particular we are interest- 
ed in a massless system of the 
type shown in Pig. 2. This 
system has the equation of 
motion : 

C6+K& - ( 3 ) 

where the system parameters C 
and K are real. 

By substitution of Eq, 2 into Eq. 3 we find: 

icjuC F + KF = F (4) 

and separation of the real and imaginary parts gives the follow- 
ing linear equations: 

-OJ F,C ^ ^K-; K 
^ I[K '- 0 

which have the solution: 

i 'z 

_ ~F / u 
6 - /C 


'\ 

(6) 



System S 



These expressions show that It is Indeed poo 
mine suitable parameters C and K If one alio 


functions of the frequency 


a - Q e 

O 



Fig. 3 

System S + m 


CO of the excltinjir for'’-^ . 

Next consider th':* dyr/i::.! ' 
system S + m, whloi. : ;• i- 

by supplementinp tj-,o 
S by a rigid maon n: j 

0 . Suppose this i'.ev; .'.vr' <,r:, 
is excited by a vertical 
harmonic force: 

Q = (y) 

acting on tho adcb tl niaor. 


Then this new cyo.f^em will i;c 

long to the class defined above and there will therefore 
exist a dlsplaoeiBent function P = + ip^ such that the 

displacement S of the mass m {and hence of the point 0 ) 
can be written in the form: 




(H 


A connection between P and F can be found by intro, h,e 
xng the reaction P acting on m (and the system S) at point 

This force is related to the displacement d through Eq. P 

and It must also satisfy the equation of motion for the 

Using the sign rule Indicated In Pig 3 tnis 
„ ^j-g. 5} this equ 

tion of motion is: 



C9) 


rr) 


I - Q 


- Re 


L Co t 


Differentiation of Eq, 2 and Eq. 8 gives: 


rn 6 


//y Cju 

k 




( 10 ) 


m 6 = - 


^ FQ, 


( 11 ) 


which, in connection with Eq. 9 yields: 

Qr, 


and 


P = 


F = 


^ k ^ 


( 12 ) 


F 


4 _ rnou^ 

K 


1 - -- r 


(13) 


The displacement of point 0 In system S + m Is conse- 
quently by Eq. 8: 


(5 = 


A y _ m LU^ 


k 


F 


(14) 


This important formula expresses the steady- state solu- 
tion to the system S + m in terms of the solution to the 
simpler system S. Its physical Importance is best seen If 
the real part is separated. This will yield the solution to 
the case when the exciting force is: 



rc'al i'i 





Performing the process of separation we find the 
placement: 



( 16 ) 


where the phase shift is: 


f ~ ~r - r"'j 

1 ^ \ 1 ■ 2 f 

The square root appearing in Eq. l6 will be cablofi 
nification factor and will be denoted hereafter b.v 
M. It is equal to the absolute value of V. 


(rr) 


t.he 

1 };<•' r:v'm!'ol 


Simple Damped Oscillator 

The use of the theory developed above is clearly rilurj- 



Fig. 4 
System S 


trated by the following simple 
example. Consider the inar.;?- 
less damped system S shown in 
Pig. 4. This system clearly 
belongs to the class of dyna- 
mic systems covered by the 
theory, and it will therefore 
have a steady-state solution cf 




for Simple Damped Oscillator 





the form: 


S-^Fe 


The equation of motion for point 0 i; 


c 8 ^ k S ^ R e 


which after substitution of Eq. l8 yields the fol]ov;!i'.r d ; s 
placement function: 


r - 




where 


IS a dimensionless number which will be called the frequnnc 
for the damped oscillator. 

Separation of the real and Imaginary parte of Kq. :o 
gives : 


1 ^ 7 / 


1 -t 


{ 22 ) 
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These two functions are shown In Pig. 5 and play an Import- 
ant role In the subsequent sections. Note that substitution 
of Eq. 22 Into Eq. 6 gives the expected result: 

C = c = constant and K = k = constant, and that the displace- 
ment function F has the asymptote: 




_ 2 
a 

o 


r 

' z 


OL. 


for a 




(23) 


Having determined the displacement function for the sys- 
tem S It Is now a simple matter to determine the displacement 
of the more general oscillator shown In Fig. 6. 

Before writing down the solu- 
tion, it Is convenient to in- 
troduce the so called mass 
ratio B defined by: 



^ t~n 


(-24) 


Fig. 6 

System S + m 


This parameter Is a dimension- 
less number which can be 
thought of as a scaled measure 
for the mass m. 


Noting that Ba^ = m oj /k we obtain from Eqs. 16, 17 
and 22 the following displacement for the general damped 
oscillator: 

(9, 


<5 = 


k 


M cos (p) 


(25) 




where 


tan (p ^ —j 


Ba 


( 26 ) 


M = 


J-Bal) +a 


( 27 ) 


The magnification factor M la of special Interest. Its var- 
iation with for different values of B Is shown In Pip. 7 
It is easily shown analytically that whenever B >1 

the response spectra shown in Pig. 7 will exhibit r. 

fa- i wxxi exhibit resonance 

peaks of the heights: 


max 


f1 = 


at the frequency ratios: 






Por B < I no peaks exist and the largest dlsni. 
for the case of static loading. 

csciii!L"i" 

found In equivalent to the more usual form 

-- •“ »- 

The advantage of using ^ = uj .. a ^ 2 

^ o ^ c/k and H = km/c^ ig ^hat 
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Fig. 7 Steady-State 
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the effect of frequency and mass can be studied separately, 
and for the systems which will be considered below this ad- 
vantage greatly outweighs the disadvantage that this method 
fails for the undamped case (c = O) . 


Elastic Half Space With Uniform Loading 

In a subsequent section we 
will need the steady- state so- 
lution for the system shown In 
Pig. 8, It consists of a per- 
fectly elastic half space 
(soil) with i;i)n mass density 
^ and the elastic constants 
G and ^ . The half space 

is excited by a vertical, 

uniform surface loading 
Lcot 



p = p^e 


per unit area, 


Pig. 8 

Uniformly Loaded 
Elastic Half Space 


and we are Interested In the 
displacement ^ of th(‘ Ih’oe 
surface. 


It is clear that it is sufficient to look at a column of 
unit area. Since no horizontal displacements occur, this 
column will behave like a rod with zero lateral displacement 
and its elastic modulus (constrained modulus) is therefore: 


Ep = G/s^ 


( 30 ) 
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where 


2 i- 2ju. 


( 31 ) 


The velocity of the P-waves propagated along the rod 
Is : 



where V is the velocity of shear waves (S - waves). No re- 

o 

fleeted wave will occur since the rod is infinitely long and 
the displacements of the rod must therefore he of the form: 

X = ^ ( 33 ) 


Which is the general expression for a sinusoidal wave with 
amplitude A propagated downwards with the constant velocity 

Vp. 

The tensile stress on a horizontal plane is: 


^(z) 




= 


oj L(ju 



( 34 ) 


and the stress boundary condition at the surface (z = O) 
gives consequently: 

/_\ • A r~ ^ 

^ - Po^ = -i e 


( 35 ) 
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or 


= 




isp^ 




fo 


uj~\^ G 


( 36 ) 


Hence, by Eq. 33, the displacement S at the surface Is: 




_ 


La/ \ 



(37) 


This is exactly the steady- state solution of the different- 
ial equation: 


1 


9^ 


6 = A ^ 


LQJ 


t 


(38) 


which is also the equation of 
oscillator shown in Pig. 9. 



Pig. 9 
Analog for 
Half Space 

displacement of the surface. 


motion for the :5!mp]e damped 

This oscillator, which con- 
sists of Just a dashpot with 
the damping coefficient: 



5 


can therefore be used a:? an 
analog for the uniformly 
loaded half space as long as 
we are only interested in the 


This analogy between energy dissipation due to wave 
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propagation and that due to viscous damping is of great in- 
terest, since it allows us to study the behavior of an infi- 
nite system by considering a simple finite system. It is 
one of the aims of this paper to show that it is possible to 
extend this idea to more complicated systems. In particular, 
we will study a system consisting of a rigid circular foot- 
ing resting on the surface of an elastic half space. 

The material elements of an elastic half space possess 
no damping properties. The observed energy dissipation is 
therefore entirely due to the infinite geometry of the sys- 
tem, and it will in the following be referred to as geometri - 
cal damping , as opposed to viscous damping , which is usually 
related to a material property, such as the viscosity of a 
fluid. 



Pig. 10 

Footing -Toil System 
It was shown on page 7 that 


We now turn to the key matter 
of the present work, the ideal- 
ized half-space model for the 
footing-soil system. This sys- 
tem, which is shown in Fig. 10, 
consists of a circular rigid 
footing of mass m which rests 
upon a homogeneous, isotropic 
and perfectly elastic half 
space with the physical con- 
stants G, ^ and ^ . 

e can find the steady-state 
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solution for this system by first finding the displace:;:-:,';. t. 
function for the case m = 0. The solution for ni 0 vrill 
then follow from Eq. l4. 

If Poisson's ratio is given, then F can at -f !r- :r, tg: de- 
pend on the parameters Cu , r^, G and g . The only dimen- 
sionless ratio which can be formed by these i s tlv- f r-'- nir-::r:y 
ratio : 


a 


UJ i 


,1 


1 




■'•01 


and since P is dimensionless we therefore conciude h.-P ;i 
is a function of and a^ only. 

The physical Importance of a^ becomes cli'ar 11’ u!;e u . 
that the velocity of shear waves (S - waves) .in an ••.las,', ic 
solid is: 



Hence the frequency ratio can also be written !n th" form;', : 


a 


0 


U) 


J? 

4 


'// 



*>, / 


Where is the wave length of S-waves in the half space-. 

These relations show that one can look at a either as 
a frequency, based on the time r /V which it takes a shear 

U Q 

wave to travel from the center to the edge of the footing, 
or as a measure for the radius of the footing in relat.ion 
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to the length of the waves propagated in the elastic half 
space . 

The static spring constant for a rigid footing Is known 
to be : 


k - 


4G 


a 


/-/V 


(^3) 


and the expression for the displacement Is by Eq. 2: 




(44) 


The parameters of the corresponding sprlng-dashpot analog 
are then by Eq. 6 and Eq. 42: 


K = 


/T 


r+ r 

'1 U 


k ^ k^k 


r = r Ks. 

' '/s ' ^ 4. 

where the coefficients c^ and k^ defined by: 

k = 


and 


— ^ 




1-2 r -2 

F +r 

'i '2 


(45) 


(46) 


are dimensionless functions of a^ and JU . 

By substitution of Eqs. 45 Into Eq. 3 we obtain the 
following equation of motion for a massless footing-soil 
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system: 




k IS 


Vs 



k,k 6 



7 ■ 


Before considering this equation in more detail we must 
study the behavior of and k^. This requires knowledge of 
the displacement function F. The analytical determi nati on 
of this function involves the solution of the wave ou/.ia'- i on 
for an elastic solid with a mixed boundary condition, name- 
ly, a zero stress condition on the free surface of the' half 
space and a uniform displacement condition under th" foot- in/’- 
This difficult mathematical problem has not yd, I”'*;; sojvMci 
exactly and it is therefore necessary to rely on ce-rta In 
analytical and numerical approximations which will he dis- 
cussed later. 


Low Frequencies 

Sung (2) and later Bycroft (3) made the plausible assump- 
tion that if the footing is small compared to the wavr‘lengl,h 

of the vibrations set up in the half space, that is if a is 

o 

small, then the pressure distribution under the foot!;;/’- is, 
close to the pressure distribution for the static case 
{ UJ = = 0). This distribution is known from the theory 

of elasticity and can be expressed by: 


Pir) 


R 


ZwrA 




( 48 ) 
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Sung ( 2 ) and Bycroft (3) therefore replaced the mixed bound- 
ary condition with the pure stress condition: 




Icjut 


= i 


2Trr^^r/-r^ 

0 


; ^ < C 


, ^ > C 


(^ 9 ) 


and were then able to obtain approximate solutions for the 
surface displacements of the half space for the low frequency 
range, aQ< 1 . 5 . 

The vertical displacement of a surface point at the dis- 
tance r from the axis of symmetry can be written in the form: 


^ ( 50 ) 

Where the distance ratio a is defined by: 

which is simply a generalization of the definition of the 
frequency ratio a^, see Eqs. 40 and 42. 

The surface displacement function g = + Igg is dimen- 

sionless and time-independent, and is characteristic for the 
chosen boundary condition. In the case of Sung's and 
Bycroft's solutions this condition is the static pressure 
distribution defined by Eq. 49. 

Sung ( 2 ) determined only the displacement at the center 
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of the footing (a = 0), while Bycroft (3) went a step further 
and calculated the displacement of the footing as a wt-^ghl.ed 
average of the displacements'" over the loaded area. Hence, 
Bycroft's solution is to be preferred over Sung's. 

Bycroft calculated the displacement of the footing posi- 
tive when upwards and wrote his final solution in i.ne foriii. 


6 



o y 


lout 


: 2 ) 


The real and imaginary parts of his dlsplacenioni- lur.ct,.)’-. 
f = f^ + if 2 are shown in Pig. 11. 

Comparison between Eq. 44 and Eq. 5H shows tha' "'ycr 
displacement function is connected to the norma 11 -d dl em- 
placement function F used in the present work througii 1 ho 
simple expression: 


f j f . 

and this formula has been used to prepare P.1g. 11' wlilrh 
shows that for any practical purpose we can con s. Id/v s' tv) h'j 
Independent of Poisson's ratio . This is, of course, the 
main reason for introducing the normalised d I spj aju 'in' cH func- 
tion. 

Prom Fig. 12 and Eq, 46, we can directly calculate the 
corresponding values of the functions c^ and . The result 
is shown in Fig. 13. 

* These displacements are not 'uniform since the asusinied 
pressure distribution is not the actual one. 
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Frequency Ratio Oo 


Fig. 11 Bycroft's Displacement Function 
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Fig* 12 The Displacement Function F 
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Fig. 13 The Functions and 
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High Frequencies 

When the radius of the footing Is much larger than the 
wave length of the waves In the half space, that Is when a^ 
Is much greater than unity, we expect the displacement to ap- 
proach the solution for uniformly loaded half space. This 
problem was solved on page l4, where we found the displace- 
ment given by Eq. 37, which can also be written In the form: 


X - ^ 

' K 


9^ 


-212 I 0~2/j) 


LOji 


(54) 


where k = 4Gr^/(l- ) = pseudo spring constant 

= an arbitrary, but large radius 
= toual force on area with radius r 


p = TTr ^ p 

O 0^0 

= (ic-' r /V 


= frequency ratio 


Eq. 54 Is In the standard form (Eq. 2) and we conclude that 
the asymptote for the displacement function P Is: 


F 


- L 


a c. 


for 


a. — ^ •=><> 

o 


where 




(55) 


(56) 


This asymptote Is similar to the one found for the simple 
damped oscillator, see Eq. 23. 
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The variation of the coefficient 1/c^ with Poisson *s 





Pig. 14 

Varation of l/coo 
with Poisson's Ratio 


ratio is shown in Pig. l4. It 
will be seen that in the range 
0.4 the error made by- 
setting 1 /cj^ = 0.92 for all 
is smaller than 6 ^ which 
is quite insignificant for most 
practical problems. That 
1 /Cfio = 0 fovyO^ = -g- is not 
surprising since this case cor- 
responds to a situation where 
the half space is incompress- 
ible. 


Intermediate Prequencies 

At very low and very high frequencies we could estimate 
the distribution of the pressure under the footing. This 
method fails completely in the approximate range 1 . 5 <'a^<^ 8 , 
where the wavelengths of the waves propagated through the 
elastic half space and the dimension of the footing are of 
the same order of magnitude. A special method will therefore 
be developed for solving the problem in this range. The 
method is based on Relssner's (l) original solution for a un- 
iform pressure distribution, and it will, for reasons which 
become apparent later, be called the ring method . 

* 

Relssner showed that the surface displacement function 


* Por definition see Eq. 50 on page 21. 
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for the case of a uniform pressure distribution over a cir- 


cular area of radius r^ is: 


/ Nx^- 


? c 
X -5 


g(a,a) = ^ 


Where 


f(x) = (2x^ - if - 4x^ix^-s^ ix^-l 
X = The positive root of f(x) 


(58) 


J & J^= Bessel functions of the first kind 
0 1 

Equation 57 is valid for all values of a and a^, but 
the Integral has never been evaluated except for a = 0 and 


a = a , for which cases it is possible to establish a solu- 
o 

tion by series. As will be understood from the following, it 
is desirable to evaluate g(a, a^) for a larger range of a and 
and such an evaluation has been carried out in Appendix 
I, for the special case of Poisson's ratio equal to 1/3. 

The actual pressure distribution under a rigid footing 
is not uniform, but varies in both phase and magnitude with 
the distance r from the center of the footing. However, in 
view of the rotational symmetry it can be assumed that the 


pressure is uniform over the area of any narrow ring which 
is concentric with the footing. The footing can therefore 
be replaced by the ring system shown in Pig. 15 . 

This system consists of n rings of equal width and the 
radii are consequently related through: 


r 


no 


nn 

n 


r 

n 


(59) 
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Fig. 15 Ring System 
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where Is Identical to the parameter r^ used above . 

The pressure over the ring is uniform and harmonic 
and can be written as: 


Pm 


R 


jr rJ 


n 


m 






( 60 ) 


where the pressure coefficient is the ratio between the 
. , th . , , , 

pressure over the m ring and the average pressure over the 
footing. 

For small amplitudes of vibration the principle of super- 
position will be valid, and the vertical displacement of 
point j, due to the loading on the m^'^ ring, can therefore be 
calculated by Eq. 50. After a series of elementary calcula- 
tions one arrives at the expression: 


where 



SQrr. 

Gr^n 





I cut 

e 


( 61 ) 


9j,. = 


(62) 


a 


m 


- 

n 


tn jr. 


CJO 


( 63 ) 


The uotal displacement of point j can be found by summation 

over all the rings, and is most conveniently written In the 
form; 



Re 


LUj 


Gen 



c„ 


'm jm 


rr):zi 


( 64 ) 
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where 


Is a known Influence coefficient . 

If the above ring system Is to be compatible with a 
rigid footing. It must be required that all the points J=l, 
2...,n have the same displacement. The following must 
therefore be satisfied: 

, J -- -l,?,-", ('"'b (66) 

Substitution of Eq. 64 into this shows that the pressure 
coefficients q^ must satisfy the n-1 equations: 

C.'- - (67) 

It must also be required that the total force over the 
n rings is equal to the exciting force . This con- 

dition yields the linear equation: 

a {Zm-1) = ( 68 ) 

J m 

which is easily obtained from Eq. 60. 

The relationships Eq. 67 and Eq. 68 form a complete set 

of n linear equations from which the pressure coefficients 

q can be determined. Both the coefficients and the 
^m 



Ay 
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unknowns in these equations are complex variables, but this 
causes little extra difficulty as the complex equations can 
be transformed Into a set of 2n linear equations between 
real quantities by simple separation of the real and imagi- 
nary parts. 

Once the pressure coefficients are determined the vert- 
ical displacement of the footing can be calculated from 
Eq. 64 with J = n, l.e.: 



_R_^ 

Gr^h 



(69) 


This equation is of the same form as Eq. 52 , and it is readily 
seen, that the displacement function f is: 

n 

i ( 70 ) 

The displacement function F follows then from Eq. 53 : 


rK,/j) 


- 4 - 


n 



m=l 


( 71 ) 


The ring method described above has been programmed for the 
IBM 7090 computer, and several runs have been made using sys- 
tems of up to 20 rings and a value of I/3 for Poisson's 
ratio. 

The most significant result is the displacement function 
P shown With a dashed line in Pig. 16 and Pig. 12 and in tab- 
ular form on page 111 of Appendix II. 


It agrees well with 
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Fig. 16 Displacement Function for Footing-Soil System 
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Bycroft's solution in the low frequency range and approaches 
the asymptote given hy Eq. 55- 

It should he noted that th e displacement function for 
the rigid footing-soil system is very similar to the dis - 
placement function for a sprlng-dashpot system , see Pig. 5. 

It was shown previously that P -is virtually Independent 

of Poisson's ratio for both low and very high values of a . 

o 

It is therefore realistic to assume that the variation of P 
is approximately as showp in Fig. I6 for all values of 
Poisson's ratio . 

Knowing P we can again calculate c^ and k^. The results 
of these calculations are shown in Fig. 17 and also^, for the 
low frequency range, in Pig. 13 (dotted curve). 

It can be shown that^^ lli ^ which was defined by 
Eq. 56. Also, it can be shown that llm k, exists, but the 

aQ— »■ oo i 

actual value is not known. This, however, is of little impor- 
tance since the term k^k S in Eq. 47 becomes Insignificant at 
high frequencies. 

Pressure Distribution 

It is not the intention here to describe in detail the 
variation of the dynamic pressure distribution under a rigid 
footing, but, since it is generated as a by-product of the 
ring method, a single example is shown in Pig, 18. The steps 
shown are the real and imaginary parts of the pressure co- 
efficients q^, defined by Eq. 60, while the curves marked 
q^ and qg are the supposed envelopes for the limiting case 
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Fig. 17 The Functions k-j^ and 
for Poisson's Ratio - l/3 




Value of Pressure Coefficients q, and qg 
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Fig. 18 Pressure Distribution for the Case: 
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of infinitely many rings. The pressure distribution is at 
any given time: 


p(rj) = 



(72) 


It Will be seen that even at the relatively low fre- 


quency ratio a^ - 1 the distribution deviates significantly 
from the static distribution, which according to Eq. 49 Is: 


7 

% 



(7.3) 


At higher frequencies the distribution becomes more com 
plicated. The central part of the distribution remains rela 
tlvely uniform while the part near the edge starts to oscil- 
late violently forming a diverging wave. 

Displacements of the Free Surface 

Another by-product of the ring method is information 
about the vertical displacements of the free surface in the 
vicinity of an oscillating rigid footing. This displacement 
will vary in both phase and magnitude with the distance from 
the footing forming the vertical component of the Rayleigh 
wave (R-wave) created by the disturbance. 

The displacement is readily calculated from Eq. 64 
which is also valid for j greater than n. This equation can 
therefore, be used to calculate not only P but also the com- 
plete surface displacement function defined by Eq. 50. Fig. 
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19 shows the real and imaginary parts of g for the special 
case, a^ = 1, and Poisson's ratio equal to 1/3- 

The displacement has the form of a divergent wave. The 
wave length, in terms of the dimensionless distance a. Is 
slightly shorter than 2 9T . By the definition of a, Eq. 51, 
this indicates that the actual wave length is slightly short- 
er than the wave length of the shear waves in the half space. 
This is in agreement with the current theories for the pro- 
pagation of R-waves. See Rayleigh (9) or for a more recent 
presentation Ewing, Jardetsky and Press (lO). 


Response Spectra 

Having determined the displacement function B’ we can 
now turn to the case where the footing has a certain mass m. 
It is here convenient to Introduce the mass ratio: 


b = 

which is such that: 



i -JU ^ 

4 ^ rj 


(74) 


Bal = mcuYK (75) 

The above mass ratio is proportional to the mass ratio 
b = m/ ^ ro used by other authors: Richart (5), Sung (2 ), etc. 

The equation of motion can be obtained from Eq. 4? simply 
by adding an inertia term: 

ml ^ c, 4^1, + 


(76) 
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which has a solution of the form: 


b 


^ /r 


(rr) 


This gives upon substitution into Eq, 7*6 the displacement 
function: 



1 

(krBa) + i-c,a^ 


US) 


The magnification factor, defined on page 8, is there- 
fore: 







(Y'O 


from which the response spectra shown by solid lines in Pig. 
20 have been calculated using the functions and shown 
in Fig. l6. 

It was mentioned on page 34 that P is practically inde- 
pendent of Poisson's ratio and the same is therefore to be 
expected for P. In fact, it is easily shown that the relative 
variation of F with is smaller than that of P. Hence, 
the spec tra shown in Fig. 20 can be used for all values of 
Poisson's ratio . 

The phase shift is: 


9 



~c,a^ 

-3a^ 


( 80 ) 



Half - space theory (/i= 1/3) 
Simplified analog (k, =1, C| = 0.85) 
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Fig. 20 Steady-State Spectra for Footing-Soil System 
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and the displacement caused by the exciting force Q^cosc<_; t 
is finally: 

S ^ M cos {cut ^ (p) (_8i ) 


Simplified Analog 

The practical application of the above theory has shown 
that only limited agreement can be achieved between observed 
and calculated amplitudes. This is mainly due to discrepan- 
cies between the theoretical half -space model and an actual 
footing-soil system, the soil phase of which is generally 
nonlinear, inhomogeneous and Imperfectly described. The 
mathematical difficulties Involved in the strict use of the 
half-space model are therefore hardly justified, and even 
crude approximations can be introduced in practical calcula- 
tions, without reducing the reliability of the calculated 
response. In particular, we can attempt to replace the half 
space system with a simple damped oscillator with similar dy 
namic properties. One way of doing this would be to assign 
appropriate constant values to the functions c^ and ap- 
pearing in Eq. 76. This can be done as follows: 

Since it has been established that c^ and k^ are bound- 
ed, see page 34 , we conclude from Eq. 78 that for B ^ 0 , 
which is generally the case, we have: 



j Vk, » 1 
1 

O 


7 


F 


1 


a, « 1 


( 82 ) 
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The last of these, with Eqs.75 and 77jgives the displacement: 


6lo i _ - Qo 

k 5al mcu^ 


(83) 


which ^hows that for B / 0 and a^» 1 the response of the 
footing-soil system is completely Independent of the proper- 
ties of the elastic half space, hence we can choose c^ and 
very freely in the high frequency range. This could also 
have been deduced from Eq. 76 by observing that the inertia 
term dominates the left side at high frequencies. 

Equation 82. shows that in order to have agreement for the 
static case we must choose = 1. It also shows that at 
very low frequencies we can assign any value to c^. 

The Influence of c^ is strongest at the resonance peaks, 
that is in the range 0.3 < < 0.8, and it is therefore 

reasonable to choose c^ as an average value over this range. 
Pig. 13 shows that O .85 is appropriate. Thus, let: 


c^ = 0.85 and k^ = 1.00 (84) 

This is equivalent to replacing the footing- soil system by 
a simplified analog consisting of a simple damped oscillator 
with the physical parameters: 
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The error made by this apparently bold approximation is 
surprisingly small as can be seen in Fig. 20 where the res- 
ponse spectra for the simplified analog are shown with dotted 
lines , 


The relative error can be calculated from: 




(k r Bcc f+[ y 
(1 -E)aff+\o.ssaS 



and Fig. 21 shows how this error varies with a^ and B. It 
will be seen that the error is small in the Important range 

0,8 where both amplitudes and resonance frequencies can 
be calculated from the analog with, an error which does not 
exceed Qfo. 

For larger values of a^ and small values of B the rela- 
tive error might be as great as However^ B-values be- 

low which correspond to very large and very light footings^ 
are not usually encountered in practice and an error of 35^ 
would therefore be most unusual. Furthermore, the absolute 
value of an error of this type will be small since it occurs 
on the tail of the response spectrum (see Pig. 20). 

Another way of comparing the solution for the half- space 
model with that of the simplified analog is to compare the 
displacement functions for the two systems. This has been 
done in Pig. l6, which corresponds to the most unfavorable 
case B = 0. 
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Practical Formulas 

Assuming the acceptance of the above simplified analog 


we can now set up the following set of formulas for the 
steady-state motion of the footing-soil system in Fig. 10: 


Spring constant 



(87) 


D _ 'A 

Mass ratio [j ~ ~7^ 

V- 


i - u. W 


( 88 ) 


Frequency ratio OL^ - 


Displacement amplitude 


Phase shift 


oj r 

o 




^o/ k 


k(i-Bcilf+{as5ay 


r - i-Bar 


(89) 

(90) 

(9.1) 


Displacement 


<S = A CCS (p) (92) 


Resonance occurs only when 0.36 and one can establish 
the following approximate formulas for the resonance condi- 
tion: 


Resonance frequency qj 


Resonance amplitude 


/! 


tes ~ 


res 


^ \B'~0.S6 


Oo 


B 

5 


(9.3) 


(9^) 


Equations 3 and 4 were arrived from Eqs. 28 and 29 by 
observing that the mass and frequency ratios for the analog 


are: 
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5^5 k,/c^ « b/o.ss" (95) 

0.65 rx^ ( 96 ) 

Free Vibrations 

Since the steady-state response spectra for the footing- 
soil system and its simplified analog show good agreement 
for all frequencies we can conclude that the two systems will 
behave similarly not only for steady-state excitation, but 
also for other types of dynamic loading. In particular, 
their free motions will be similar. 

The natural damped frequency for a simple linear oscil- 
lator is: 




k S ~ 0.25 
^ . 


(97) 


Hence, by Eqs. 85 and 95 j we obtain the following expression 
for the natural damped frequency of a footing-soil system: 


w,x = 


ii? -Q.is 
S 


(98) 


which is good for B ^ l/'^* B-values smaller than this 

the motion is either non-osclllatory or very strongly damped. 
Such small B-values are rarely encountered in technical appli- 
cations . 

The logarithmic decrement for a damped oscillator can 
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be shown to be: 





~0.zs 


(99) 


where and ^ ^ ^re consecutive displacement amplitudes 
to the same sl^g. gy substitution of Eq. 95 Into this we 
obtain the foliov^ing approximation to the logarithmic decre 
ment for a footing-soil system: 


^Ln ^ - O.id ( 100 ) 

Mass ratios higher than B = 5 are rarely encountered in 
practice and si^ce the logarithmic decrement for this case Is 
as high as 1 . 22 , corresponding to x^ / x^^^ = 3.4, we con- 


clude that Qillj:^rmal fo oting-soil systems are strongly damp - 
ed by geometrical dispersion of wave energy. 

y important dynamic systems, such as mechanical os- 



Pig. 22 

Rotatlng-Mass System 


dilators and machine founda- 
tions, are excited by one or 
more eccentric rotating masses. 
A typical member of this class 
of systems is shown in Pig. 22. 
The system is identical to the 
system S + m in Pig. 3 except 
fop the mass m^ which rotates 
relative to the mass m with 
sccentricity ^ and angular 
frequency 60 . The spring 





50 


constant K and the damping constant C are frequency dependent 
and are defined by Eq. 6 . 

The equation of motion for the mass m Is: 

(m+mj 6 + C6 A 6 m^coswt (loi) 


which shows that the system behaves like a system with the 
mass m + m^ excited by the harmonic force cos co t. 

For the special case of the simple damped oscillator , 
see Pig. we can write this force as: 

Q ^ cos cut ( 102 ) 

where 

" (4j 


Substitution of Eq, 102 Into Eq. 25 gives directly the 
displacement: 

= -f- K M co^{cut + (p) (10_4) 

where M and (p are given by Eqs. 26 and 27 respectively. 
The variation of the magnification factor: 



— 2 
a. 





where 

B = (106) 

is Shown in Pig. 23 . por B > J It exhibits a resonance 
peak at:=, •; ,, ■» 



(107) 
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with 

/nax M, = iXfd - 'A ( 108 ) 

and for all values of B it approaches the limit l/B as the 
frequency goes to infinity. This implies that the amplitude 
of displacement has the limit: 


/ Im 6 


J nif 
ni + 


(109) 


cu oo ' ' ' ' ■ 

This last formula is true for all rotating -mass systems 
of the general type shown in Fig. 22. 

For the case of the half-space model it is practical to 
Introduce the mass ratio: 


O _ fT) + nrij ( Vs \ _ / -A> Wiotal 

^ K V^j 4 ^ rj' 


(110) 


With this variable the equivalents of Eqs.102 to 105 become: 


a - 

<5 = 




cos Lui 


1 


( Vs. V 




alM = 




(in) 

( 112 ) 

(113) 

(114) 


The phase shift cp is, as for the constant force oscillator, 
given by Eq. 80. 

wp.fV# 
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The variation of is shown in Fig. 24. The full 
lines correspond to the solution for the half-space model 
(c^ and from Fig. 1?) and the dotted lines correspond to 
a new simplified analog obtained by setting: 

k^ = c^ = 0.90 (115) 


for all frequencies. This choice of k^ and c^ gives a some- 
what smaller error in the practical region of resonance, see 
below, than the simplified analog defined by Eq. 84. 

The new simplified analog is a simple damped oscillator 
for which: 


“ Vi//^ = ( 116 ) 

B = 5 k^cf ~ B/o.s ( 117 ) 

These relations, in connection with Eqs. 107 and 108 give 
directly the following approximate formulas for the deter- 
mination of the resonance peaks: 


a = 
0 


A 


o.s 


5 -0.45 


maK 




0,9 


B -0.225 


(118) 

(119) 


No peaks exist when B ;5 0.45 and the maxima for 
B O .9 are hardly significant enough to be classified as 
resonance peaks. Even for the heaviest footings we rarely 
encounter mass ratios greater than 5, and the practical 
range of resonance frequencies is therefore approximately 
0.7 < < 1.4. 
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Fig. 24 Rotat ing-Mass Spectra for Footing-Soil System 



Ill TRANSIENT MOTION 

Numerical Method 

The transient motion of a strongly damped linear system 
excited by a oulse of arbitrary shape and duration can be 
determined by a simple Fourier analysis if the steady- sta’te 
solution for the system is known at all frequencies. In 
order to show this we first consider the motion of a footing- 
soil system excited by the vertical dynamic loading shown in 
Pig. 25. This loading consists of a train of identical 
pulses of alternating sign and duration d. The pulses are 
equally spaced with the time interval Jd (j = integer), and 
the loading is therefore periodic with the period 2jd. The 
individual .pulses can be expressed in the form: 

a(i) = 

where is the maximum value of the force, and ^ is a 
dimensionless pulse-shape function , see Fig. 26a. 

The pulse train shown in Fig. 25 can be developed into 
a Fourier series of the form: 

N 

Qit) - ( 1 . 1 ) 

/?=0 

where N is a large integer and 
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Period = 2jd 


I 


H 


Fig. 25 Typical Pulse Train 


0(t/d) 



Fig. 26a Fig. 26b 


Typical Pulse-Shape Functions 
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( 1 - 22 ) 


X - 

cx„ = ^ j <p(u) cos ^ du 

y.„ = fjjiu) .Ln(^u)da ' 


(123) 

(124) 


For small amplitudes the principle of superposition is 
valid and a Fourier series for the displacement can therefore 
be obtained by applying Eq. 8 l to each of the terms of Eq. 
121. This gives: 


where 




a 





8n-^f 

tan^ 



K~ 


(125) 

( 126 ) 

(127) 

( 128 ) 


The variables and c^ are functions of a^ and can be 
obtained from Fig. 17 or calculated from Eq. 46. Or, if it 
is desired to base the calculations on the simplified analog 

they can be set at the constant values 1.00 and O .85 respec- 
tively. 
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Equation 126 shows that: 



Bal 



const, ^ for 


t) 


x.: ( 129 ) 


hence the series defined by Eq. 125 converges faster than 
the series given by Eq. 121. This means that a good approx- 
imation to the response of the footing can be obtained with 
a reasonably small number N of, harmonics of the pulse train. 
This corresponds to the well-known phenomena for multi-degree 
of freedom systems^ that most of the response occurs in the 
lower modes of vibration. A method for choosing the number N 
will be given below, see Eq. l43. 

Now consider the behavior of the system under the in- 
fluence of a pulse train with, large spacing between the 
pulses (j large). During the period while a pulse is "kcting, 
energy is fed into the system and the amplitude of vibration 
will therefore increase. During the subsequent period, when 
no force is acting, the system will perform free vibrations 
which, since the system is strongly damped, will die out 
after a. few oscillations. Hence, when the consecutive pulse 
arrives, the system will behave as if the previous pulse had 
never existed. 

Consequently, if the above analysis is performed with a 
sufficiently large choice of j, the result will be a periodic 
displacement which, within each half period corresponds to 
the displacement due to a single pulse. 

The above calculations are easily performed on an elec- 
tronic computer and a FORTRAN II program for pulse shapes of 



the piecewise linear type shown in Pig. 26 Td is presented in 
Appendix II. The program ^calculates not only the displace- 
ments but also the velocities and the accelerations which are 
easily obtained by differentiation of Eq. 125. 

The piecewise linear pulse used In the program is partic- 
ularly convenient because it makes it possible to replace the 
Integral expressions of Eq. 124 by the following finite sums 
which are easy to handle in the computer: 



yi-, y-L*) 

i-~o ^ 

{zn^fTT^ H ^^yry-i-r y^i) 

^ i~o 


cos 

sin 




2n^i 

jm 


IW 


r (130) 


where 


y-,^ y^^ym^r ^ 

The pulse type also has the advantage that the exciting 
force can be defined simply by specifying the ordinates y^, 
^2"‘" ^m-1" duration d and the force level Q^. 

The Integers N and j 

A few remarks should be given as to the choice of the 
Integers j and N. 

If the footing-soli system is of the oscillatory type^ 
B large ^ it has^ according to Eqs. 98 and 100, the approxi- 
mate natural period: 



^rr 




(132) 
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and the logarithmic decrement 


A 


C’ 


in 




(133) 


Hence, the amplitude of free vibration will vary as; 

- (2.7 t \ 

/i = vi^r/ 


(134) 


and the time t^ necessary to reduce the amplitude to a 
fraction 8^ of the initial amplitude is: 

r 

Loi 

's 

For footing-soil systems of the non-oscillatory type, B 
small, it can be shown by considering the simplified analog 
that the amplitude varies according to the law: 


S' ^7 ' ■■ !/ 


“ -SB-^logS, (135) 


A = A e 


/JJi) 


(136) 


hence, we have for this case; 


/ = - dS5 -jj- in 

"i i/ ^ 


8-^ Log 
' 10 


137) 


The estimates, Eq. 135 and Eq. 137^ can be combined into the 
upper bound: 


t ~ -{8+ 6 B)^ Log f 
^ Vs io ^ 


(138) 


which for all values of B can be used as an appropriate time 
interval between the individual pulses of the pulse train. 
This gives the following criterion for the integer j: 


> y - (2 3- 6S; ^ 


J 


(139) 



bO 


An appropriate choice of the Integer N can be obtained 
by the condition that the magnification factor corresponding 
to the highest frequency ratio a^^ appearing In the series Eq 
125 should be smaller than a certain value^ say . This 

gives by Eqs, 79 and 84 the following restriction on a^^: 


L / c2 


(1^0) 



The Integer N can often be chosen considerably smaller 
than Indicated by this Inequality but such a reduction requires 
special consideration In each case. 

It Is usually rational to choose the ^3 S-ppearlng In 
Eqs. 139 and l 43 to be Identical, say 0.05 or 0.01. Smaller 
values are hardly justified for practical calculations. 



TrapezoldaX- Pulse 

Consider the footing-soil system defined by: 


K, ~ ft , G — 2500 psL 
IV = 403 kips , X = 128.&B pcf h 

. ylC = ^ 

The velocity of shear -waves in the half space is by Eq, 4l: 

1 / i 2 500 ^ i44 iy± / ^ . f 1 ) \ cz\ 


1/ = \ - ^nn Pf /<pr (i^5) 

\| f 26 . 6 B/s 2 J 7 ri/set 

and the mass ratio is by Eq. 74: 


B = 


(/ -h)^ 432>,000 
iZS.ib X 5^ 


(146) 


The latter shows that the system is of the oscillatory t3^< 


The damped natural frequency is by Eq. 98- 

200 ~A5~0m _ o/ _ -/ 

ou , - P - Zb.?> sec 

oL o o 


(147) 


and the corresponding period is: 


T = -^2 ^ 0.24 sec 


( 148 ) 


The logarithmic decrement is by Eq. 100; 
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i.ZZ 


-aid 


(149) 


and the ratio between consecutive amplitude maxima in free 
vibration is therefore 3.4 indicating a strongly damped 


system. 

Now let this system be excited by the pulse shown in 
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Pig, 27. This pulse has the property of being constant, 
equal to 250,000 lb., for a period of O.9 sec., which is 
considerably longer than the natural period of the system. 



Pig. 27 Trapezoidal Pulse 


Drawing from our knowledge of linear damped oscilla- 
tions we therefore expect the footing to perform damped 
oscillations about the static displacement: 

a, ^ ^ " 

4-^ 2500 ■< 144 ^ 5 ^ 



^sictiLc 


As a check on this expectation the response was calcu- 
lated by the method outlined above. Thp 6 's were chosen 
to be 0.05 which gave: j = 2 and N = 40. The actual compu- 
ter output Is shown on page 116 , and the displacements are 
shown graphically In Pig. 28 (full curve marked B = 5, 

G = 2500). It will be seen that the natural frequency, the 
amount of damping, and the static displacement agree with 
the figures calculated above. 


The curves marked with the mass ratios B = 2, 1 and ^ 
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show the effect of changing the weight of the footing while 
all other variables are being kept constant. The free 
motion of a light footing is, as expected, more strongly 
damped than that of a heavy footing. 

An Increase of the modulus G of the soil to say 
G = 5000 psi affects the magnitude of the vibrations of the 
system, but not its relative damping. This is Illustrated 
by the curve marked B = 5^ G = 5000. 

All the above curves were obtained by applying the 
above Fourier method to the half -space model; l.e., using 
the functions c^ and k^ shown in Fig. 17 in the formulas 
Eqs. 126 and 128. 

Exactly the same method can be used for the simplified 
analog defined by Eqs. 85. This is achieved simply by set- 
ting c^ = 0.85 and k^ = 1 and the result for the case B = 5, 

G = 2500 psi is shown with a dotted line in Fig. 28. The 
close agreement between the response of the simplified analog 
and that of the footing- soil system shows, as expected, that 
the analog can be used, not only for the case of steady-state 
motion, but also for transient motion. 

This observation is of importance for practical calcula- 
tions since the analog is a classical system, which can be 
treated by a variety of well known methods, such as numerical 
integration or phase-plane techniques, which do not require 
the availability of an electronic computer. 
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Triangular Pulse 

As a second example, we 



Pig. 29 

Triangular Pulse 


will consider the behavior of 
a footing-soil system excited 
by the simple triangular pulse 
shown in Pig. 29 . The dis- 
placement is given by Eq. 125, 
which after introduction of 
the dimensionless time: 



and the corresponding dimensionless duration: 



can be written in the form: 


6(r) = 


(152) 


(153) 


This expression shows that the maximum displacement can be 
expressed by: 

'it 

Where D ( t 1 B) is a dimensionless function which will be 
called the displacement spectrum . Its variation with ^ 
and B is shown in Pig. 30, which was determined by application 
of the above Pourier method to a series of selected footing- 
soil systems excited by triangular pulses of different dura- 
tion. The dotted curves in the figure show corresponding 
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values calculated from the simplified analog. The close agree- 
ment between the two set of curves confirms again the feasi- 
bility of the suggested analog. 

Displacement spectra of the type shown in Pig. 30 can 
be prepared for any pulse shape, not only for maximum dis- 
placements but also for maximum velocities and accelerations. 
They are of great practical interest since they constitute a 
simple and fast tool for the calculation of the maximum re- 
sponse of a given footing-soil system. 



IV COMPARISON WITH TEST RESULTS 


Several researchers have presented data which substan- 
tiate the elastic half-space theory presented above^ at 
least for small displacements. Most of the data was obtain' 
ed from steady- state tests in the low frequency range, 

< 1.5, and special attention was given to the resonance 
condition. Lately this data has been supplemented by re- 
sults from Impact- loading tests. 

Determination of Soil Properties 

In studying the test results one must bear in mind 
that only limited agreement with the theory can be expected 
since soil is not a perfect elastic material. And, even 
under laboratory conditions It Is not homogeneous nor Iso- 
tropic. The determination of appropriate elastic constants 
to be used In the theory Is therefore a difficult problem. 

An extensive survey of the available methods has been given 
by Whitman (ll) and only two of the more promising methods 
will be described here. Both methods are dynamic, which Is 
natural since the dynamic shear modulus of soil can vary 

widely from the static value, and both methods involve the 
determination of wave velocities. 
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The first method Is trlaxial dynamic testing; . This 

technique has been developed by Shannon et al. (12)^ Hardin^ 

Rlchart (13) and Hall (l4). It involves the determination 

of the torsional and longitudinal resonance frequencies of 

cylindrical trlaxial specimens. Prom these frequencies one 

can calculate the wave velocities V and V and the elastic 

s p 

constants of the half space follow then from: 


and 


n = 
^ 



(155) 


( 156 ) 


Equation 155 can be derived from Eqs. 31 and 32 and Eq. I 56 
follows directly from Eq. 4l. 

The tests must be carried out at a confining pressure 
corresponding to the average pressure under the footing. 

This can be estimated by the theory of elasticity. For 
clean sand one can obtain a first estimate of G from the 
test results published by Hardin and Richart (13)- 

The second method to be discussed here will be called 
in-situ testing . It consists of two parts. First, the vel- 
ocity of P-waves is measured by a simple seismic test. In 
this test a shock is induced at the soil surface by a drop 
weight or hammer and the arrival time of the first distur- 
bance is recorded by a pick-up placed on the soil surface at 
a known distance from the source. Since the P-waves travel 
considerably faster than any other type of waves this gives 
a direct measurement of Y^. By comparing the arrival times 
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at different distances It Is furthermore possible to deter- 
mine the variation of with the depth below the surface. 

For further details see Pry (17). 

The second part of the In-sltu method Is a determination 
of the velocity Vj, of Rayleigh surface waves by a method sug- 
gested by Jones (15) and further developed by Ballard (16) 
and Pry (17). in this method a small vertical oscillator Is 
placed on the soli surface and the nodes of the resulting 
surface waves, R-waves, are located by means of a portable 
pick-up. Prom the distance between the nodes and the known 
frequency one can then calculate the velocity of the R-waves. 
This velocity Is only slightly smaller than the velocity of 
S-waves. For soils one can assume that = 1.05 7 . The 
actual ratio varies with Poisson's ratio and Is equLalent 
to the variable defined on page 28. Knowing we can 
again calculate Poisson's ratio ^ from Eq. 155 and the 
Shear modulus Q from Eq. 156. Without going Into details It 
Should be mentioned that the method can be refined to detect 
variations In with the depth below the surface. This Is 
done by performing the test at different frequencies. The 

Idea being that low frequency R-waves penetrate deeper than 
waves of high frequency. 

Steady- State Tests 

Steady-state tests on footings were carried out long 
before the appearance of the first elastic half-space theor- 
ies. Notable are the tests carried out at the Deutsche 
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tests were analysed by Hertwlg, Pruh and Lorenz (l8) who^ 
like the author, compared the vibrations of the footing- soil 
systems with that of a simple oscillator. However, while the 
author's model has the same mass as the footing, the above 
researchers were trying to set up a set of rules for deter- 
mihlng the so called in-phase mass , which is the mass which 
has to be added to the mass of the footing in order to obtain 
agreement between the observed and the theoretical response. 
The search was only partly successful since it turned out that 
the' in-phase mass depended on frequency and several other 
factors. In spite of this, the concept of the in-phase mass 
has since achieved great popularity among engineers and 
Hsieh (4) has shown how a frequency-dependent in-phase mass 
can, in fact, be calculated from the half-space theory. 

The U. S. Army Engineer Waterways Experiment Station 
at Vicksburg, Mississippi, has carried out a very extensive 
program of steady-state field tests on circular concrete 
footings having diameters varying from 5 to l6 feet and ex- 
cited by an oscillator of the rotating-mass type. Data were 
obtained for 63 tests on a silty clay (loess) and 39 tests 
on a uniform fine sand. The results were reported by Pry (6) 
and later analysed by Chae (7) and Whitman (ll), who both 
found good agreement between the observed resonance ampli- 
tudes and the corresponding amplitudes calculated by the 
elastic half- space theory. The standard deviation from the 
predicted amplitudes was about 30^ which must be character- 
ized as satisfactory under field conditions. As could be 
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expected the agreement with the elastic theory was best for 
the cohesive soil, 

A series of steady-state tests of special Interest to 
the author was carried out at the University of Michigan by 
Chae (7) whOj simultaneously with the author's theoretical 
work, determined the displacement functions f^ and f^ exper- 
imentally and studied the dynamic pressure distribution 
under a rigid footing on dry sand. For this purpose, a 
special footing was constructed, see Pig. 31. it consisted 
of a set of concentric steel rings, each of which was con- 
nected to a common rigid steel plate through four load cells. 


The latter were simply short pieces of aluminum tubing in- 


strumented with SR-4 strain gages. These, In connection 
with additional electronic equipment, made it possible to 
determine the pressure distribution under the lower plate 
due to a static or dynamic loading on the upper plate. 

The soil used In the tests was highly compacted Ottawa 
sand contained In a 4 '-9" square sand bln with. 4' high con- 
crete block walls. The system was driven by an oscillator 
which could be placed on the upper plate of the footing by 
means of a special loading frame. The mass ratio B of the 
system could be varied from approximately 1 to 5 by means of 


special weights. 

A sugary of Chae's test results for the displacement 
function P is shorn, In Pigs, 32 and 33. He actually found , 
small variation with the mass ratio of the footing, hut thl, 
has been neglected for the purpose of the present 


presentation . 
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AREAS: 

RING I TO 5 * 21 in’ 
ring 6 « 8.1 IN* 



Fig. 31 Special Footing, from Chae (7) 
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The agreement with the author's theory Is good except for 
at very low frequencies. This is not surprising since it 
is well known that the static pressure distribution under a 
footing on sand deviates drastically from the theoretical 
distribution given by Eq. 49. The actual pressure has a 
maximum under the center of the footing and vanishes at the 
edge. The effect of this change on the elastic half- space 
theory was studied by Sung (2) who assumed that the pressure 
distribution was parabolic and obtained the solution for P^ 
shown in Pig. 32. It should be noted that Sung's solution 
is for the center of the loaded area (r = 0). The average 
displacement is therefore somewhat smaller than indicated 
by the curve. The large deviations on Pg found at high fre- 
quencies can be explained by experimental errors^ in partic- 
ular difficulties in determining the phase shift and the fact 
that at high frequencies the Inertia force from the mass of 
the footing over shadows the forces from the soil. These 
latter are therefore difficult to determine accurately and 
this is again reflected in the large scatter of the experi- 
mental values for P^ and Pg, The inertia effect is, of 
course, the same as the one from which we profited on page 
43 when we established the simplified analog. 

Impact Tests 

Data on Impact tests is scarce in the literature but a 
series of tests which support the author's theory has been 
carried out at the University of Michigan and the results 
will be published shortly by Drnevlch, Hall and Rlchart (8). 



Parabolic distribution (Sung) 
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.xperimental Values lor , by Chae 
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Fig. 33 Experimental Values for F- , by Chae (7) 
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The footing and the soil were exactly the same as that used 
by Chae (7 ) but the oscillator was replaced by a mechanism 
with which a small sand bag could be dropped on the top of 
a load cell placed on the upper plate of the footing. This 
simple arrangement made it possible to record the actual im- 
pact load on the footing and experience showed that both the 
shape and the duration of the pulse could be controlled by 
changing the height of fall and by placing different foam 
rubber pads on top of the load cell. The displacements were 
recorded in terms of accelerations which were obtained by 
means of an accelerometer mounted on the footing. The re- 
sults of two tests are shown in Pigs. 3^ and 35 where they 
are compared with the theoretical response calculated by 
the method described in Chapter III. The actual computer 
Output for one of the curves in Pig. 35 ( =1/'^) is shown 
on page 123 of Appendix II. 

The agreement between the theoretical and the experi- 
mental curves is excellent up to a certain time, approximate- 
ly 37 nis for test M-1, and then the curves suddenly lose 
their similarity. This can be explained as follows: The 
principal effect of the Impact is the propagation of an 
S-wave from the footing towards the side walls of the sand 
bln and a P-wave towards the bottom of the bln. Both of 
these waves will be reflected by the walls. The half-space 
theory is therefore valid only up to the time when any one 
of these reflected waves arrive back at the footing. This 
return time is easily estimated. Por test M-1 the shear 



wave velocity is = 325 ft/seo and assuming Poisson's ratio 
equal to 1/4 we find, by Eq. 32 , = 560 ft/seo. The S-wave 

has to travel approximately 4.75 ft; thus it returns back in 

15 ms. The P-wave has to travel 8 ft and it returns to the 

footing in 14 ms. Since the peak of the pulse occurs at ap- 
proximately 18 ms we can therefore expect serious deviation 
from the theory after 33 ms. This number agrees well with 
the 37 ms actually observed. For test T-3 the waves travel 
faster, therefore the disagreement occurs somewhat earlier. 

Finally, it should be noted that it makes little dlff- 

erence which value of Poisson's -pai-no -? =, ,<1 • 

roisson s ratio is used In the calcula- 

tions. Therefore, it is hardly worthwhile to carry out 
apeclal tests to establish a correct value for this ratio. 

An intelligent guess should be sufficient for most practical 
calculations. A value of 0.33 for dry soil and 0.50 for sat- 
urated soil will usually be accurate enough. 



Footing: ro= 0.5 ft, W = 69.7 lb 
Soil: G = 2500psi, X = 109 pcf 
Observed 

Calculated, /i= 0.25, B = 0.96 
Calculated , fj, - 0.40 , B = 0.77 


Fig. 34 Test M-1, by Drnevich, Hall 


and Richart (g) 
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tn 
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Footing: ro=0.5ft, W = 239.9 lb 
Soil: G* 4200 psi,X= 109 pcf 
Mass Ratio: 8*3.30 


0 ms 


Observed 

Calculated , 0.25 



4 20 


\\/ 


time 


Fig. 35 Test T-3 
by Drnevich, Hall and Richart (fi) 



V CONCLUSION 


The half-space model considered In the preceding 
chapters is a compromise between mathematical convenience 
and accurate representation of an actual footing-soil system. 
However, the usefulness of the model has been clearly demon- 
strated by the fact that the discrepancies between calculated 
and' observed displacements fall well within the range which 
one can expect for a soil-footing system, which by nature 
is rather poorly defined. 

On the other hand, the observed discrepancies are sig- 
nificantly larger than the differences found between dis- 
placements calculated from the half-space model and corres- 
ponding displacements calculated from the simplified analog. 

The conclusion to be drawn from this must be that the simpli- 
fied analog is to be preferred for its mathematical ease over 
the -only apparently- more accurate half -space model. Thus, 
the problem of determining the vertical motion of a footing - 
soil system has been reduced to the problem of determining 
the motion of a simple damped oscillator defined by the follow- 
ing equation of motion: 


mS ^ <5 ^ = Qit) ( 157 ) 
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This analog, which is the most important .result of the pre- 
sent work, can be solved by several well known methods such 
as direct or numerical integration, phase plane, etc. 

The study shows that all footing-soil systems are 
strongly damped due to wave propagation into the subsoil. 

This geometrical damping will generally be much more signifi- 
cant than the material damping of the soil and the latter can 
therefore be neglected in practical calculations. The damp- 
ing ratio decreases with Increasing mass ratio of the system 
and the motion of a large and light footing will therefore 
be more strongly damped than that of a small and heavy foot- 
ing . 

The theory developed in the preceding chapters is a 
linear theory. It therefore applies only as long as the dis- 
placements can be considered small as compared to the size 
of the footing. Also, it is assumed that the footing and the 
soil are in physical contact at all times. This means that 
the soil reaction R(t) = P(t) + W must not become negative in 
which case the footing will Jump free of the ground. This 
latter limitation is of little consequence for practical de- 
sign, since footings are usually dimensioned with a certain 
safety factor against Jumping. However, several steady- state 
test results. have been published which Indicate that Jumping 
occured during the tests and in some of these tests the devia- 
tions from the half -space theory were considerable. 

Only the vertical mode of vibration has been considered. 
However, the nature of the energy dissipation for the other 
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modes (rocking, sliding, torsion) is basically the same as 
that of the vertical mode, and it is therefore reasonable to 
expect that simple spring-dashpot analogs can be developed 
also for these modes. Once these analogs have been developed 
it will be a relatively simple matter to include the dynamic 
characteristics of the foundation in the analysis of a super- 
structure with several degrees of freedom. 



appendix I 


REISSNER'S DISPLACEMENT FUNCTION 
Relssner's surface dlsplacenient function was defined 
on page 28 but will be repeated here for convenience: 

oo 


Where 


= ( 1 - 2 ^ )/2(i-/y ) 
/< = Poisson's ratio 
f(x) = (2x^-1 


2 2 
X -s 


1 


2 

X -1 


x^ The positive real root of f(x) 
and , Bessel functions of the first kind 
The function f (.) la called H ^ylelgh function It 

roo"t ; -- positive 

Hthouf oPiW. 

ithout going into details it should h 

mentioned that x 

veioc rr Tr --- - 

The independent variables a and a ar<- 
their defini^- ' according to 

nition, always positive and real a 

integration parameter x. ^ 
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The integral appearing in Pn icp • 

q. 1^8 IS not directly suited 

for numerical integration for the foil • 

^oe following reasons: 

1) The integrand Is not defined In the point x = x 

Where f(x) has a root ° 

2) The interval of Integration is Infinite. 

3 ) The integrand is in general complex. 

") % = a = 0 the integral can attain 

infinitely many values, depending on the manner 
in Which a and a^ approach the origin. 

Definitions 

Fortunately all the above die-Pnr, 

difficulties can be overcome. 

n or er to show this It Is convenient first to consider the 

behavior of the following functions, which are all related to 

the Rayleigh function: 


fM = 


f%) 


(Zk-J) + ^ 


0^ 

/ ^ X < OO 

0 ^ K ^ S 
f ^ K < Co 


( 159 ) 


( 160 ) 


r (x) = z'w/'fx) = xU(2s^-3) x^- 0^./ i 161 ) 


The intervals shown with Eqs. 159 and 160 are the In 
tervals In which the functions attain real values. 



Also consider the following functions: 


r{x) = 


>6es(x) j 0 ^ K ^ S (163) 


t(x) = jbesM i ( 164 ) 


u(x) = 


(2x^-/) /'('x) jbesM x<oo (165) 


ix6<) = -jx f (x) jbes(x) , / ^ X < 00 (166) 


U/(x)= ix(x) + ^(^X^-x) ^ / X x<oo (167) 

yP^x) = /* Cx^/fX-Xo) j 0 ^ X < CO (168) 


h(x) =*1x^-5^ f(x) bes(x-l/P(x) = 1/ (x) (x-x^^) (169) 


The Intervals shown with the definitions are the inter- 
vals in which the functions attain real values. The question 
as to whether the functions are actually defined in all 

points of the given Intervals will be left to the subsequent 
sections. 

It Is clear that the function f*(x) has no root greater 
than s. Hence, the function f»»(x) has only one root x = x 
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in the Interval x > s. Eq. l68 shows that, since f**(x) is a 

. . 2 p p 

polynomial in x and has the root x = x . also P(x) is a 

o 

2 

polynomial in x , Furthermore , since P(x) has no roots for x 
greater than s the function h(x) is finite for all values of 
X greater than unity. 

By substitution of Eq. I69 into Eq. 16? we obtain: 


or 


\J!/ (X) 


h (k) h (ZXq-x) 


(170) 


IV/ (X) 


(x-3Xo)/?fx) + (K + Xo)h(2 Xo-x) 

(x-xj(x + xj(x-3xj 


(171) 


This function is not defined in the point x = x^ where 
both the numerator and the denominator have a root, but its 
limit for x-*-x^ exists and can be found by L 'Hospital's 
theorem. In order to show this we differentiate the numera- 
tor and the denominator of Eq. 17I in the point x = x^: 


num'i>^o) ~ Z h(X-o) “ 

4 x^h'(x 4 ) 

(172) 

den' (K o) - ~ 4 x^ 


(173) 


and L 'Hospital's theorem gives directly: 

lim M/(X) = h(x,) - (174) 

Since h(x) is finite and has a finite derivative in the 
point X = x^ it is obvious that the limit exists. Hence, the 
singularity of w(x) in the point x = x^ is removable and w(x) 
is continuous for all values of x greater than unity. The 
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actual value of the limit can be determined, but it is of 
no interest as the calculations can be planned in such a way 
that it is not necessary to know this quantity. 

For the special case s = -g both f*(x) and f**(x) have 
2 

the root x = and the functions r(x), t(x), and u(x) are 
therefore not defined in this point. However, also these 
singularities can be removed, which is easily shown by 
L 'Hospital's theorem. 

From What has been said up to this point it is a simple 
matter to establish, that the following functions are real, 
finite and continuous in the Intervals stated below: 


/•(x) 

r 

0 

4 

X 

4 

S 

t (x) 

1 

0 

4 

X 

4 

/ 

^(x) 


5 

4 

X 

4 

i 


2Xo 

-1 

4 

X 

< 

CO 


) 

1 

4 

X 

4 



( 175 ) 


Separation of the Real and Imaginary Parts 

Relssner ' s displacement function can be written as: 


Sf ^ 9, * % ( 176 ) 

where g^ and gg are the real and Imaginary parts respectively 
By simple complex algebra, which will not be repeated here, 
it can be shown that: 

oo 

V(x)afK ( 177 ) 




89 


and 


9 

'-'Z 


^jr(x)dx ^ ± 


tM c/k 


" be;.(x) ( 178 ) 




where, according to Eq. 175^ all the integrals are real. 

The Integrand v(x) is not defined in the point x = x^, 
but this difficulty can be overcome by writing: 


. oo 


"2 A.-/ 


OC> 


V (x) dK - I l/'(x}dx + I \/(k)o/x + / '/ dx 

4x,-/ 

■Xo r ^ 

= I [^^(x) - v^(2xo-x)] c/x + / k (x) dx 

4x^-; 

- I W(x)dx + / v(x)dK + TT p( 


(179) 




where 


R 



(180) 


is the remainder after cutting off the Integration at some 
large value K of the parameter x. 

Substitution of Eq. 179 into Eq. 177 gives: 




Xo 




= ^L(x)dx + ^\w(x)dx-^ ^ i/(x)c/x + /? 


(181) 




In this formula, and in Eq. 178, all the integrands are 
real and continuous in closed intervals of integration and 
the Integrals can therefore be evaluated by numerical integra- 


tion if X and K are given, 
o 
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Calculation of x and f'fx ) 

Q i 

The root x^ Is most conveniently determined from f**(x) 
by Newton's Iterative method: 

xL- (182) 

where 

M ^ 4d(s‘-l)x -16(2s^-3)x‘-d ( 183 ) 

This will yield the values of both x^ and d(f**(x^) )/d(x^), 
and f'(x^) can then be calculated as follows: 

By Eqs. 159 and I6O: 

fk)- 1(2x1- if . 3 ,, 


By differentiation of Eq. I6I: 


d 




(185) 


and finally by Eqs. l 84 and I85: 



d . 

cJM' ( 186 ) 


.•^.^y^pbot 1 c Expression for v(x) 

In order to determine the "out-off" K and the remainder 
R it is necessary to study the behavior of v(x) for large 
values of x. According to the definition of v(k) we have: 
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VM = bes(x) (187) 

which is finite and real for x > x^'^ 1 . It can he shown by 
standard methods that: 

_ o _ o 

where o(x'’'^) stands for a term which decreases as x~^ when 
X and a similar expression can be obtained for the 
function bes(x) = J (ax) • J, (a x) if the Bessel functions 
are replaced by their asymptotic expansions: 




cos 




+ 


sin jcLx-'^) 
da X 


(189) 

-f 




Wax 

o 


cos (a,x- - 


3 sin (cLoK - W) 


Sa^X 


+ o(x^-^) 


(190) 


which are valid for large values of ax and a x. 

o 

By multiplication of Eqs. I89 and 19O we find: 

COS (cLX~%) OOS {a^K 


bes(x) - 




+ 


sin (ax-%) cosjcLoX—^^^^ 


Sax 


_ 3 cos (ax ~% ) sin 


da^X 


(191) 
+ o (x-^) 
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which can be reduced to: 


bes (x) 


TTX-iaZ. 


sLn(cL^-a)x - cos(o^a'^a}K 

+ (ac~a)x - sin (a^-^a)K 

&a.K (1$ 

_ cos (ao-a)x + sZr? (fgp+gjy + q 


Substitution of Eqs. l88 and 192 into Eq. I87 gives 
finally the following asymptotic expression for v(x): 


l/(x)' 


where 


^ Zin AfCL]x - C(n:> (a.^ + a)x 

+ cos(a^-a)x -f- sln(a.^-i-a) K 


(193) 


27T 


( 194 ) 


c = ^ 

" daa 


'195). 


3a - qLq „ 


(196) 


THe, expansion Is not valid for a or a^= 0 , and these cases 
will therefore have to be Investigated separately: 

For a ^ 0, but a^^_£_0, the above method can be used to 
show that the asymptotic expression for v(x) is: 


Kx) « tS 




(197) 
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where 






(198) 


c 


-3 


da. 


(199) 


The case a = a =0, which corresponds to the case of 

Q 

static loading, will be discussed on page 102. 


Determination of K and R 

Substitution of the above asymptotic expressions into 
Eq. 180 gives the following formulas for the remainder R: 
For a 7^ 0, a^ ^ 0, 1: 

'OO fOO 


R 


J_ 

TT 


CA )R-sln(cL„-a)xclx - cJA^cos(a/a)KdK 


JK 

oo -^>00 


+ 


cj x^cos(a^~a)KdK -t c jx^sin(a^-\a)xdK 


'K 




(200) 


For a = 0, a^ ^ 0, K 1; 

.00 


R ~ ^jc/x + -^jK'%in(a,X-^dx {201) 

'^K 

Both these expressions contain integrals of the types: 


- lx^cos(xx~lAu)dK 


(202) 
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and 



((XX - ihir)dK 


(203) 


where m is greater than unity. These integrals can be eval- 
uated as follows: 

First consider the case CQ 0 . The integral can then 
be written as: 



oo 

(o<x-/b<!r)dK 


(204) 


where 



= I x'^cos((XK -/b9r) dx 


(i^+fh)ir/c 


0< 


and is chosen in such a way that: 


(205) 


(206) 

Where is a large Integer. 

The substitution x = u + (n + /». )V-/k reduces Eq. 205 
to : 




1 

II 

u + (n^^)7r 


( 207 ) 
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since n + ^ > 0 this shows that the terms of the series in 
Eq. 204 are of alternating sign. By a well known theorem it 
can therefore be concluded that the absolute value of the 
series is smaller than its first term. This value can be es- 
timated if we rewrite Eq. 207 in the form; 




n m-f 



h-^^7T)cosu - (u + slnu 


du (208) 


from which it can be seen that: 


7 / (A?y-/s) - (n^/iArj) . m (X (209) 




Hence, by Eq. 206: 




11 

< 



< 


m IT 


1 


( 210 ) 


This result can be generalized to cover also the case 
when (X < o, provided Eq. 206 is replaced by: 


r 


+ 4 ) <K^ (211) 


1(X| 'b( 


•/ “ |o(| 'oc 


Equations 204 and 210 show that the error made by assuming 
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that : 


-K, 

1 = I x'”cos(c<K-ii 7 !')c/x , <x 4 0 


is numerically smaller than: 


( 212 ) 


A 


/ 


rn'll 

2‘“T7m4T 


D( 


K 


(213) 


provided is chosen according to Eq. 211, 

Similarly it can be shown that the error made by re- 
placing Ig with: 



)('^6ln(o(.K 


fbTT^olK 


iR 


(/.¥ 0 


is numerically smaller than: 


(214) 


^2 


(215) 


provided Kg is chosen such that; 


+ (ai6) 


where Ng is a large integer. 
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Finally, by noting from Eqs. 211 and 2l6 that: 

K-K < (- 7 ) 

we conclude from Eq. 212 that the error made by neglecting 
or Ig entirely is smaller than: 

OT ~rr) 

A = t7 A' (278) 

-'3 |e<| 

Now let and Kg be two large numbers, which are 
chosen in such a way that: 



TT 





irNi 
« +a 

o 


(219) 


where is taken as the smallest integer which satisfies 
this condition. Then Eq. 200 can be written in the form: 


■K, 


f\ + -^ \ ){cos{a^-^a)'Kdx — 


(220) 


7r 


- 00 


r 00 

c I x^sin (cX- ot}xc/K - cos(a^djxd>( 


2 

00 


' c>o 

f-3. 


+ ^^cos[a-d)KdK + cJk un[a^a)xd){ 

L "A' K, 
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Hence, by Eqs. 213^ 215 2l8: 


R + ^ x^cos (a +a)xdx 

k 


< 


+ 


2* |c.| , 2\c\ 

{ci-afK\ ia,^a)K] 

jCal ICal 

\a-a\K\ 


( 221 ) 


which, in connection with the approximation ~ Kg, shows 
that the error made by putting: 


rK. 


R - 


V 

^ I v72 

7r 


x^cos(a+a))(dK 




is smaller than: 


(222) 


A,= a; 


-3 

/ 

a 



2Q 



1 

CTs 

4_ 

C3 



Ldfaf 

T- 

{a+a)^ 


a-a 


oL^-^a 



(223) 


or by Eqs. 138-140: 


+ 


A. = 


{a^~af (a^+a) 


3cl O-o 

daa^(a„-a) 


+ 


3a -a. 


daa^(a^) 


2r(/-s‘)-^^x: 


(224) 


It is a simple matter to verify that A|j_ is smaller 

than: 




( 225 ) 
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where Is the smallest of the quantities a, ' 

And this expression for the error will be used in the follow 

Ing as it is much more convenient to work with than Eq. 224. 

It is clear from Eq. 224 that the approximation Eq, 222 

cannot be used for the calculation of R when a = a , A spec 

o ^ 

ial formula is therefore needed for the case a = a ; 

o 

Substitution of a = a^ into Eq. 200 gives: 


oo 


oo 


oo 


.A' 


-^Ix^Cos2clkc/x 4 - 


^jx dx+ -p X%lo 2cLKdx 


( 226 ) 


whichj after evaluation of the second Integral, reduces to: 




R 


Cz 




kS 


OO 


^ jx\os 2axdx + ^oLxdx i^^T) 


Now if is chosen according to Eq. 219j. it can be 
seen from Eqs. 213 and 2l8 that the error made by assuming 
that : 


R 


2W^i 




( 228 ) 


is numerically smaller than: 



-3 


K. 


' iQl 




2cl 


(229) 
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or by Eqs. 19 ^ and I 96 : 


^ 


(230) 


which is similar to Eq, 225. 

The remainder for the case a = Oj a^ 7 ^ 0 is given 
by Eq. 201 which in connection with Eqs. 215 and 2l8 shows 
that : 


A 


7 



sIqI 

o 


+ 



or by Eqs. I 98 and 199 : 


(231) 


A 


7 


iS 



0.7S 


(232) 


provided the "cut-off" is chosen such that: 



M + 


TT 


(233) 


where is a large integer. 

Formulas for the Calculation of g ^ and 

Prom the above expressions for the remainder it is 
now a simple matter to set up procedures from which the fun- 
ction g^ can be calculated to any given accuracy A . Depend- 
ing on the values of a and a^ we must distinguish between the 
following three cases: 
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1) a^O, / 0, a 7^ a^: 

Choose two large K- values such that: 


y K ^ JT Ni _ y IZ ^ TJA 


where 


a^ = The smallest of the quantities a,a^j |a-a^| 

II2 = A large integer 

= The smallest integer which satisfies Eq. 234, 


and calculate from: 


r! r 

I (/(x)c/x f / ^x(:x)dK + 

— iwlZrlx + / co s,(cl + a, )xc/x 


^ = 


L £a'.*-/ 

2) a^ ^ 0, a 0, a = a^: 


Choose a large K- value such that: 


(235) 


\ < K ^ 

d%S0-s>)A 


fa36) 


where N is a large Integer, and calculate g, from; 


f r° 

UMc/K+ W(k)g/k + 


1 i/(x)c/x - 


avHi-d) cdK^ 


’2k.-1 


(237) 
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3) a = 0, ^ 0: 

Choose a large K-value such that: 


1.Z 




of <K = 


where N is a large integer, and calculate g^ from: 

f /Xo rK 

+ -|r /l^(x)tyx + -~jv(X)c/K 




(238) 


( 239 ) 


5 '/ 

The function can for all values of a and be calcu- 
lated directly from Eq. I78, which will be repeated here for 
convenience: 










fO^o) 


bee>(x) ( 240 ) 


All the functions appearing in the above formulas are de- 
fined on page 85 and they are well behaved in the intervals 
of interest. The integrals can therefore be evaluated by nu- 
merical integration. 


The Static Case 

Reissner's displacement function is not defined in the 
point a = a^ = 0, which according to the definition of a 
and a^ corresponds to the case of static loading. 

The vertical displacement of the surface for a uniform 
static loading of intensity p^ over a circular area is known 



to te : 






—P c 

/ ti o 


TT (1- s") C- ■ 1 


jVz 

cos\ cIk 

- sin^K 

■% 


^ > C, 

5 / o 


(241) 


y-(^)'rfx, /•< 






By comparing these with the following form of Eq, 50: 




(242) 


it can he seen that if CU-^0, while the ratio 

kept constant, the components of g(a,a ) will 

o 

limits: — , 


Cos^x dx 


L irn 0 
O' -^(9 


- 1 

7r^(f-sV 


I -^(cL/af-sw^x 

Vz 


^0 


'y_(^)V, 


a/a^ = r/r is 
' o ' o 

approach the 

^ ay etc 

(243) 

) CL <a^ 


I Im a = 0 ( 244) 

CC' Q '-^2 

The elliptical Integrals appearing in Eq. 243 cannot he re- 
duced further, hut it is known that their values lie in the 

range 0 to 1 depending on the ratio a/a^. Hence, the function 

2 2 

g^ can attain any value in the range 0 to l/TT (s -l) in the 
point a = a^ = 0, depending on the ratio between a and a^ as 
they approach zero. 

The limits Eq. 243 are useful if it is wanted to tabu- 
late g^ for small values of a and a^. Instead of using a 



104 


and a as entries In the table one can use a and a/a if 
o o 

a < a or a , and a /a, if a > a . By this procedure Is 
0 O' i 

transformed into a well behaved, single valued function and 
the resulting tables are well suited for interpolation. 

Numerical Calculations 

A computer program was written for the procedures out- 
lined on page 100 and the functions g^ and g^ were evaluated 
in the range O^a^^ 8 and 0^a^8 for Poisson's ratio equal 
to 1/3, using the formulas Eqs. 178 to l84. A graphical 
presentation of the results Is given in Pigs. 36 and 37. 

Also a table of the type mentioned in the previous sec- 
tion was prepared on punched cards. These cards made the 
values of Reissner's function easily available, and they 
formed part of the input for the computer program for the 
ring method. See page 32. 
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ace Displacement Function 
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CO 



CM 

a» 


The Imaginary Part of Reissner’ 
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APPENDIX II 

COMPUTER PROGRAM FOR IMPULSIVE LOADING 

The computer program shown on the following pages Is 
written in the FORTRAN II language. Only the most basic 
statements have been used and the program should therefore 
be able to run on most modern computers with only minor 
modifications . 

The program calculates the dynamic response of the 
footing-soil system shown in Fig. 10 due to the piecewise 
linear pulse shown in Pig. 26b. The method used is the 
Fourier technique described in Chapter III. 

The comment cards shown on page 109 explain how the 
input should be arranged and two sets of output are shown 
on pages ll6 and 123. These outputs correspond to the 
curves shown in Pigs. 28 and 35 respectively. 

It should be noted that the program automatically gen- 
erates an appropriate time Interval for the final table 
showing the actual response. This time interval will al- 
ways come out as a "nice" number for easy plotting of the 
response curves. The program also determines the number m 
of Intervals in the pulse by locating the last non-zero 
ordinate from the data. Hence, the number m is not part of 
the data as one might expect from the development in Chapter 
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III. 

The controllng variables E . and EPSl and EPS2 in 

the program, are preset to the common value O.O 5 . They can 
be chosen smaller, but experience with the program has shown 
that only a small gain in accuracy is achieved by doing so. 

As a special feature the program checks the sign of the 
soil reaction: 


R(i) = Qii) ( 2 '* 5 ) 

where g = acceleration at gravity, to establish whether jump- 
ing occurs. If R(t) is negative the acceleration in the 
final table is printed out with an E-conversion Instead of 
the usual E-conversion. 

The subroutine shown on page ll4 generates the displace- 
ment functions and Eg for the half -space model. These 
functions are shown graphically in Fig. 16 and numerically 
in the table on page 115 . The functions were evaluated for 
Poisson's ratio equal to 1/3. However, as shown on page 34, 
the same subroutine can be used for all values of Poisson's 
ratio. 

The execution time varies with the number of points nec- 
essary to define the pulse, the duration of the pulse and the 
choice of the preset variables EPSl and EPS2. The execution 
times for the two data sets presented below were approximate- 
ly 30 seconds each. 
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C FORTRAN II PROGRAM TO CALCULATE THE DYNAMIC RESPONSE OF A 
C RIGID CIRCULAR FOOTING RESTING ON AN ELASTIC HALF-SPACE 
C 

C CODED BY JOHN LYSMER* UNIVERSITY OF MICHIGAN* MARCH 1965 
C 

C THE PROGRAM CONSISTS OF A MAIN PROGRAM PLUS A SUBROUTINE 
C WHICH YIELDS THE DISPLACEMENT FUNCTION F = FI + I*F2 
C 


c 

EACH DATA 

SET CONSISTS 

OF 11 

CARDS AS FOLLOWS 

c 

CARD 1 





c 

COLS. 

1-10 

R 

( FT ) 

RADIUS OF FOOTING 

c 

COLS. 

11-20 

W 

( LB) 

WEIGHT OF FOOTING 

c 

COLS. 

21-30 

GAMMA 

(PCF) 

DENSITY OF SOIL 

c 

COLS. 

31-40 

G 

(PSD 

SHEAR MODULUS OF SOIL 

c 

COLS. 

41-50 

MU 


POISSONS RATIO OF SOIL 

c 

COLS. 

51-60 

D 

( SEC) 

DURATION OF PULSE 

c 

COLS. 

61-70 

Q 

( LB ) 

FACTOR ON PULSE ORDINATES 

c 

COLS. 

71-80 

IDENTIFICATION. MAY BE LEFT BLANK. 

c 



ANY CHARACTERS PUNCHED IN THESE COLUMNS 

c 

r 



WILL 

APPEAR 

AS A HEADING ON THE OUTPUT 

c 

CARD 2 





c 

COLS. 

1-10 

Y( 1 ) 

FIRST 

NON-ZERO ORDINATE OF PULSE 

c 

r 

COLS. 

11-20 

Y( 2 ) 



c 

COLS. 

71-80 

Y( 8 ) 



L 

c 

AND SO 

ON UP 

TO 



L 

c 

CARD 11 





c 

COLS. 

1-10 

Y{73) 



c 

r 

COLS. 

11-20 

Y(74) 



V-. 

c 

COLS. 

71-80 

Y( 80) 




C 

C THE DATA CAN BE PLACED ANYWHERE IN THE RESPECTIVE 
C FIELDS PROVIDED THE DECIMAL POINT IS PUNCHED. 

C ALL BLANK FIELDS ARE ASSUMED TO MEAN ZERO. 

C 

C THE OUTPUT CONSISTS OF 

C 1. A TABLE OF THE DISPLACEMENT FUNCTION USED 

C 2. A PAGE SHOWING THE DATA AND CONSTANTS OF THE SYSTEM 

C 3. A TABLE SHOWING THE MAGNITUDE OF THE FOURIER TERMS 

C 4. A TABLE OF DISPLACEMENTS. VELOCITIES. ACCELERATIONS 

C AND THE EXCITING FORCE (FROM THE TRUNCATED FOURIER 

C SERIES). WHENEVER THE ACCELERATIONS ARE PRINTED IN 

C THE FORMAT XXXX.X INSTEAD OF THE USUAL X.XXXE+XX IT 

C MEANS THAT THE FOOTING JUMPS FREE OF THE HALF-SPACE 

C AND THE RESULTS MUST BE TAKEN WITH SOME RESERVATION 

C 

C THE CONSTANTS EPSl AND EPS2 ARE PRESET TO 
EPSl = .05 
EPS2 = .05 

C THEY CAN BE CHANGED BY CHANGING THE TWO PRECEEDING CARDS 



no on on 
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DIMENSION Y(85)*C{ 400)»PSI( 400)»BN( 400)»PHASE( 400) 

TABLE OF DISPLACEMENT FUNCTION 
A = 0. 

CALL TABLE(-1.,F1,F2) 

1 CALL TABLE( A »F1*F2) 

print 100»A»F1*F2 

A = A+.2 
IF(A-8,1) 1,1,2 


READ INPUT FROM THE 
2 READ 
READ 


11 DATA CARDS 

114»R,W»GAMMA*G»P0IS»D*Q*A1»A2 
101»(Y{I)» 1=3.82) 


C 

C 


C 

c 


CALCULATE BASIC CONSTANTS OF THE SYSTEM 
V = SQRTF(G/GAMMA*4632.5) 

SPRING = 576.*G*R/( l.-POIS) 

RV = R/V 
RVD = RV/D 

B = W/32.17/SPRING/RV/RV 

AJ = 2.- L0GF(EPS1)*( .85+2.6*B)*RVD 

J = AJ 

AJ = J 

IF(B) 4.4*3 

3 AN = l.+AJ/6.28/SQRTF(B*EPS2)/RVD 
N = AN 

AN = N 

4 PJ = 3.14159265/AJ 

FIND THE NUMBER OF INTERVALS IN THE PULSE 
DO 6 I = 3*82 
IF(Y(I)) 5.6.5 

5 M = I 

6 CONTINUE 
Ml = M-1 
AM = Ml 
DT = D/AM 


print data and constants 

SPRING = SPRING/12. 


PRINT 
PRINT 
Yd) = 0. 

Y(2) = 0. 

T = 0. 

PRINT 

M2 = M + 1 
Y(M2) = 0. 

DO 7 I = 2. M2 
LINE = 1-2 
FORCE = Y(I)*Q 
PRINT 


113.A1.A2 
102. R.W.G. GAMMA 


108. D.Q 


.POIS.V. SPRING. B 


103.LINE.T.Y(I), FORCE 



n n n n on 
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IF(LINE-26)7»36»7 

36 PRINT 112 

7 T = T+DT 

IF (( LINE-21 )*( LINE-25) ) 42.42*43 

42 PRINT 112 

43 PRINT 111.EPS1.EPS2 

CALCULATE FOURIER CONSTANTS FOR THE EXCITING FORCE 
PRINT 109 

N1 = N+1 
DO 16 K = l.NI 
AK = K 

H = PJ*(2.*AK-1. ) 

ALFA = 0. 

BETA = 0. 

DO 8 I = l.M 
AI = I 

COEFF = 2.*Y( I+l )-Y( I )-Y( 1+2 ) 

ANGLE = H*(AI-1*)/AM 

ALFA = ALFA+COEFF*COSF(ANGLE) 

8 BETA = BETA+COEFF*SINF{ANGLE) 

C(K) = 2.*Q*AM/AJ/H/H*SQRTF(ALFA*ALFA+BETA*BETA) 
SWITCH = 1. 

CALCULATE PHASE ANGLE 

9 IF(ABSF(ALFA)-ABSF(BETA) ) 11.11.10 

10 PHI = ATANF(-BETA/ALFA) 

IF(ALFA) 12.13.13 

11 PHI = 1.5707963-ATANF{-ALFA/BETA) 

IFCEETA) 13.13.12 

12 PHI = PHI+3. 14159265 

13 IF(SWITCH) 14.15.14 

14 PSI { K) = PHI 

CALCULATE FOURIER CONSTANTS FOR DISPLACEMENTS 
A = H*RVD 

CALL TABLE{ A.F1.F2) 

FF = F1*F1+F2*F2 
ALFA = F1/FF-B*A*A 
BETA = -F2/FF 

BN(K) = C(K)/SQRTF(ALFA*ALFA+BETA*BETA)/SPRING 
SWITCH = 0. 

GO TO 9 

15 PHASE(K) - PSI(K)+PHI 
IF(K- 49). 16.41.37 

37 IF(K-103) 16.41.38 

38 IF(K-157) 16.41.39 

39 IF{K-211) 16.41.40 

40 IF{K-265) 16.41.16 

41 PRINT 112 

16 PRINT 104»K.A.C(K) .BN(K) 
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C CALCULATION OF TIME INCREMENT 
IF(B-.5) 18*18»17 

17 DT = RV/5.*SORTF(B) 

GO TO 19 

18 DT = MAXlF{D/AM/2. ♦D*AJ/100. ) 

19 ALOG = L0GF{DT)*.43A294A82+18. 

CHR = INTFCALOG) 

ALOG = ALOG-CHR 
IF(AL0G-.84) 20»20»25 

20 IF(ALOG-.50) 21.21»24 

21 IF(AL0G-.16) 22. 22*23 

22 DT = 10.**(CHR-18. ) 

GO TO 26 

23 DT = 10. **(CHR-17. 69897) 

GO TO 26 

74 DT = 10. **(CHR-17, 30103 ) 

GO TO 26 

25 DT = 10.**(CHR-17. ) 


PRINT FINAL TABLE OF DISPLACEMENTS ETC. 

26 T » 0. 

27 JUMP = 50 

PRINT 105 

28 IF(JUMP) 27,27,29 

29 JUMP = JUMP -1 
IF(T-D*AJ) 30,30,35 

30 ACCEL = 0. 

DISPL * 0. 

VLOCT = 0. 

FORCE = 0. 

DO 31 K = 1,N1 
AK = K 


31 

32 

33 

34 

35 


H = PJ* (2.*AK-1,)/D 
BNK = BN(K) 

ANGLE = H*T+PHASE(K) 

DISPL = DISPL + BNK* COSF(ANGLE) 

VLOCT = VLOCT - BNK* S I NF ( ANGLE ) *H 
ACCEL = ACCEL - BNK* COSF ( ANGLE ) *H*H 
FORCE = FORCE + C ( K ) *COSF ( H*T+PSI ( K ) ) 

IF ( FORCE+W*( 1 . -ACCEL/ 386. 4 ) >33,32,32 

1 06, T, DISPL, VLOCT, ACCEL, FORCE 

107, T, DISPL .VLOCT, ACCEL, FORCE 


110 


PRINT 
GO TO 34 
PRINT 
T = T+DT 
GO TO 28 
PRINT 
GO TO 2 


C 
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100 FORMAT (21X,F3*1,2F15.5) 

101 FORMAT (SFIO.O) 

1020FORMAT( IHl/// 14X,13HF 0 0 T I N G/ 14X » 1 3 ( IH- ) / 17X ♦ 

16HRADIUS* 19( IH, ) »3HR =F9.2*4H FT/ 17X fSHWEIGHT ♦ 19 ( IH. ) 
2.3HW =F9.2.4H LB//14X.7HS 0 I L/14X»7 ( 1H-) / 

317X.13HSHEAR MODULUS ♦ 12 ( IH. ), 3HG =F9*2*5H PSI/17X. 

428HDENSITY ...GAMMA =F9.2»5H PCF/17X* 

528HP0ISS0NS RATIO MU = F9.4/17X, 

628HSHEAR WAVE VELOC I TY. . . . . .V =F9.2»8H FT/SEC//14Xt 

711HS Y S T E M/14X*11 ( IH- )/17X, 

828HSTATIC SPRING CONSTANT. ..K = F9.0*9H LB/INCH/17X. 
910HMASS RATI0»15 ( IH. ) ,3HB = F9.4/ ) 

103 FORMAT ( I 2 1 » 2X , 3 ( 1 PEI 2 . 3 ) ) 

104 FORMAT ( 1 18 .3 ( 1PE14.3 ) ) 

1050FORMAT (IHl/// 11X.32HT I M E DISPLACEMENT VELOCITY 
122H ACCELERATION FORCE // 13X.3HSEC* 7X, 

241HINCHES INCH/SEC INCH/SEC/SEC LBS / ) 

106 FORMAT ( 7X ♦ 5 ( IPE 12 . 3 ) ) 

107 FORMAT ( 7X ♦ 3 ( 1PE12 . 3 ) . OPF 12 . 1 » 1 PE12 . 3 ) 

1080FORMAT (14X,9HP U L S E / 14X»9(1H-)/ 

117X» 28HDURATION =1PE10.3»5H SEC / 

217X. 28HFACTOR ON ORDINATES =1PE10.3»5H LBS // 

318X,41HP0INT TIME ORDINATE FORCE / 

428X.5H(SEC) ♦19X.5H(LBS) /) 

1090FORMAT (IHl/// 24X»28HF OURIER SERIES/ 

124X»28(1H- )/ 

215X.47HTERM FREQUENCY AMPLITUDE OF FOURIER TERM/ 

3I4X.48HNUMBER RATIO FORCE DISPLACEMENT 

4/17X» 1HK8X»4HA( K) ♦lOX.AHCIK) * lOX » 5HBN ( K ) / ) 

110 FORMAT (IHl / / / / lOX » 1 8HLOOK FOR MORE DATA ) 

1110F0RMAT( 1H0»13X» 15HA C C U R A C Y / 14X ♦ 15 ( IH- ) / 

117X,6HEPS1 =F6.4.33H GOVERNS DISTANCE BETWEEN PULSES/ 
217X.6HEPS2 =F6.4»33H GOVERNS NUMBER OF FOURIER TERMS) 

112 FORMAT (IHl // ) 

113 FORMAT (IHl/// 34X»2(A5)) 

114 FORMAT (7F10.0»2(A5) ) 


END 
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SUBROUTINE TABLE ( A ♦ FI , F2 ) 

IF(A) 1»2.3 

1 PRINT 100 

1000FORMAT( 1H16(/)21X*33HT ABLE OF FI AND F2 
l/S21*33(lH-)/23X» 31H HALF-SPACE MODEL / 

2 23X» 31H POISSONS RATIO = 1/3 / 

3/22X.lHA.12X»2HFl»13X*2HF2 /) 

RETURN 

2 FI = 1. 

F2 = 0. 

RETURN 

3 Y = A*A 

IF( A-1.) 4»4*5 

4 FI = l.-Y*( .60000-Y**15000) 

GO TO 6 

5 FI = (.25000 +( 1.10000-.80000/A)/A) /Y 

6 IF(A-1.18) 7,7t8 

7 F2 = -A*( .78378-Y*{ .25800-Y*. 03030 ) ) 

RETURN 

8 F2 = (-.95493+( ,7873-(l. 618-1. 276/A)/A)/A)/A 
RETURN 

END 


Note: 

The above subroutine is a polynomial approximation 
to the displacement function F determined by the ring 
method described in Chapter II. It can be used for all 
values of Poisson’s ratio. 
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T A 

B L E U F FI 

A 

NO F2 


HALF-SPACE MODEL 

POISSONS RATIO = 1/3 

^0 

FI 


F2 

.0 

l.COOOO 


.CGOCO 

.2 

.97624 


-.15470 

.4 

.90784 


-.29731 

.6 

.80344 


-.41690 

.8 

.67744 


-.50486 

1.0 

. 55000 


-.55608 

1.2 

.42438 


-.57003 

1.4 

.32018 


-.53791 

1.6 

.24414 


-.48961 

1.8 

. 18957 


-.44341 

2.0 

. 15000 


-.40314 

2.2 

. 12081 


-.36888 

2.4 

.09886 


-.33979 

2.6 

.08206 


-.31495 

2.8 

.06898 


-.29357 

3.0 

.05864 


-.27501 

3.2 

.05035 


-.25874 

3.4 

.04363 


-.24437 

3.6 

.03810 


-.23159 

3.8 

.03352 


-.22014 

4.0 

.02969 


-.20982 

4.2 

.02645 


-.20047 

4.4 

.02369 


-.19195 

4.6 

.02133 


-.18416 

4.8 

.01929 


-.17700 

5.0 

.01752 


-.17040 

5.2 

.01597 


-.16429 

5.4 

.01462 


-.15861 

5.6 

.01342 


-.15333 

5.8 

.01236 


-.14840 

6.0 

.01142 


-.14379 

6.2 

.01058 


-.13947 

6.4 

.00982 


-.13540 

6.6 

.00914 


-.1 3157 

6.8 

.00853 


-.12795 

7.0 

.00798 


-.12454 

7.2 

.00747 


-.12130 

7.4 

.00701 


-.11823 

7.6 

.00659 


-.11532 

7.8 

.00621 


-.11255 

8.0 

.00586 


-.10991 
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TRAPEZOID 


FOOTING 



RADIUS 

..R 

— 

5.00 

FT 


WEIGHT 

• « W 

— 

483000.00 

LB 

s 

0 I L 






SHEAR MODULUS 

..G 


2500.00 

PSI 


DENSITY GAMMA 

= 

128.68 

PCF 


POISSONS RATIO 

.MU 


.3333 



SHEAR WAVE VELOCITY.... 

..V 

= 

300.00 

FT/SEC 

s 

Y S T E M 






STATIC SPRING CONSTANT. 

..K 

sr 

9COOOO 

LB/INCH 


MASS RATIO 

..B 

= 

5.0047 



PULSE 


DURATION 


= 1. 

OOOE CO 

SEC 

FACTOR 

ON ORDINATES. 


= 2. 

50CE C5 

L8S 

POINT 

TIME 

ORDINATE 

FOR 

C E 


(SEC) 



(LBS) 

0 

.OOOE 00 

.0006 

CO 

.OCCE 

CO 

1 

5.000E-02 

l.OOOE 

00 

2.500E 

C5 

2 

lO.OOOE-02 

l.COOE 

00 

2.5G0E 

G5 

3 

1.500E-01 

l.COOE 

CO 

2.5C0E 

05 

4 

2.000E-01 

l.OOOE 

00 

2.5C0E 

05 

5 

2.5C0E-Q1 

l.OOOE 

00 

2.5C0E 

05 

6 

3.000E-01 

l.OOOE 

CO 

2.5C0E 

C5 

7 

3.500E-01 

l.COOE 

00 

2.5CCE 

05 

8 

4.C00E-01 

l.OOOE 

00 

2.5G0E 

C5 

9 

4.500E-01 

l.OOOE 

00 

2.5C0E 

G5 

10 

5.000E-01 

l.OOOE 

00 

2.5GCE 

05 

11 

5.5C0E-01 

l.OOOE 

00 

2.50CE 

05 

12 

6.000E-01 

l.OOOE 

00 

2.5G0E 

05 

13 

6.500E-01 

l.OOOE 

00 

2.5C0E 

05 

14 

7.000E-01 

l.OOOE 

00 

2.5CCE 

C5 

15 

7.500E-01 

l.OOOE 

00 

2.5C0E 

05 

16 

8.0COE-01 

l.OOOE 

CO 

2.500E 

05 

17 

8.500E-01 

1 . OOOE 

00 

2.500E 

05 

18 

9.000E-01 

l.COOE 

00 

2.5C0E 

05 

19 

9.500E-01 

l.COOE 

00 

2.5C0E 

05 

20 

lO.OOQE-01 

.OOOE 

00 

.OCOE 

CO 

A C C U R 

A C Y 






EPSl = .0500 GOVERNS DISTANCE BETWEEN PULSES 
EPS2 = .0500 GOVERNS NUMBER OF FOURIER TERMS 
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F 0 U R I 

6 R S E 

R I 

TEKF 

FREQUENCY 

AMPLITUDE UF 

NUMBER 

RATIO 

FORCE 

K 

AIK) 

CIK) 


I 

2.6186-02 

2.1606 

05 

2 

7.854E-02 

8.313E 

04 

3 

1.309E-01 

3.5146 

04 

4 

1.8 336-01 

3.9186 

04 

5 

2.3566-01 

1.450E 

04 

6 

2.8806-01 

2.6316 

04 

7 

3.4036-01 

6.3616 

03 

8 

3.9276-01 

1.9636 

04 

9 

4.45IE-01 

2.C41E 

03 

10 

4.9746-01 

1.5236 

04 

11 

5.498E-01 

5.299E 

02 

12 

6.0216-01 

1.195E 

04 

13 

6.5456-01 

2.1046 

03 

14 

7.0696-01 

9.3376 

03 

15 

7.5926-01 

3.0306 

03 

16 

8 . 1 166-01 

7.186E 

03 

17 

8.639E-01 

3.5C0E 

03 

18 

9.1636-01 

5.3966 

03 

19 

9.6876-01 

3 * 640 1 

03 

20 

1.021E 00 

3.9106 

03 

21 

1.0736 00 

3.538E 

03 

22 

1.126E 00 

2.6956 

03 

23 

1.1786 GO 

3.264E 

03 

24 

1.230E 00 

1.7266 

03 

25 

1.283E 00 

2.8766 

03 

26 

1.335E 00 

9.7956 

02 

27 

1.3B8E 00 

2.4236 

03 

23 

1.440E 00 

4.3476 

02 

29 

1.492E 00 

1.946E 

03 

30 

1.545E 00 

6. 7136 

01 

31 

1.597E 00 

1.478E 

03 

32 

1.649E 00 

1.4866 

02 

33 

1.702E 00 

1.0456 

03 

34 

1.754E 00 

2.395F 

02 

35 

1.806E 00 

6.6876 

C2 

36 

1.859E 00 

2.3306 

02 

37 

1.911E 00 

3. 6026 

02 

38 

1.963E 00 

1.5626 

02 

39 

2.0166 00 

1,262E 

02 

40 

2.0686 00 

3.461E 

01 


s 


FOURIER TERM 
DI SPLACEMtIMT 
BN ( K ) 

2.4C8E-C1 
9.5136-02 
4.245E-02 
5.159E-C2 
2. 163E-C2 
4. 6646-02 
1. 4166-02 
5.596E-C2 
6. 2046-03 
3. 511E-G2 
8.406E-C4 
1.357E-C2 
1.801E-C3 

6.2796- 03 
1.651E-C3 
3.2496-03 
1.339E-C3 
1.7726-03 
1.04CE-03 
9.8246-04 
7.885E-04 
5.3696-04 
5.852E-C4 
2.8C2E-C4 
4.255E-04 
1. 327E-04 
3.018E-C4 
4.9976-05 
2.C71E-C4 
6 • 640E— C6 
1.3616-04 

1. 2796- 05 
8.422E-C5 
1.810E-C5 
4. 753E-C5 
1.560E-C5 
2.276E-05 
9.328E-C6 
7. 135E-C6 
1.8566-06 



f^4 fsi m 


TIME DISPLACEMENT 


SEC 

-OOOE GO 
l.OOOE-02 
2.000E-02 
3.000E-02 
4.000E-02 
5.000E-02 
6.000E-02 
7.000E-02 
8.000E-02 
9.000E-02 
l.OOOE-01 
1. lOOE-01 
1.200E-C1 
1.300E-01 
1.4G0E-01 
i.500E-0i 
1.6COE-01 
1.700E-01 
1.800E-01 
1.900E-01 
2.000E-01 
2. lOOE-01 
2.200E-01 
2.300E-01 
2.40CE-01 
2.500E-01 
2.600E-01 
2.7C0E-01 
.800E-0 1 
.9C0E-01 
.OOOE-01 
. lOOE-01 
3.200E-01 
3.300E-01 
3.400E-01 
3.500E-01 
3.600E-01 
3.700E-01 
3.800E-01 
3.900E-01 
4.000E-01 
4.100E-01 
4.200E-01 
4.300E-01 
4.4UOE-01 
4.500E-01 
4.600E-01 
4.700E-01 
4.800E-01 
4.900E-0 I 


INCHES 

4.276E-05 
1.729E-03 
6.839E-03 
1.847E-02 
3.947E-02 
7.123E-02 
1.130E-01 
1.620E-01 
2. 143E-01 
2.664E-01 
3.142E-01 
3.55ie-01 
3.870E-01 
4.087E-01 
4.196E-01 
4. 198E-01 
4. 105E-01 
3.932E-01 
3.696E-01 
3.421E-01 
3. 128E-01 
2.838E-01 
2.57iE-01 
2.342E-01 
2. 164E-01 
2.044E-01 
1.984E-01 
1.983E-01 
2.036E-01 
2.134E-01 
2.266E-01 
2.420E-01 
2.584E-01 
2.746E-01 
2.a96E-01 
3.023E-01 
3. 122E-01 
3. 189E-01 
3.222E-01 
3.222E-01 
3. 192E-01 
3.137E-01 
3.063E-01 
2.977E-01 
2.a85E-01 
2.794E-01 
2.711E-01 
2.640E-01 
2.584E-01 
2.547E-01 
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VELOCITY ACCELERATION FORCE 

INCH/SEC INCH/SEC/SEC LB 


8.345E-02 
2.924E-01 
7.831E-0 I 
1.592E 00 
2.632E 00 
3.707E 00 
4.599E 00 
5.144E CO 
5.281E 00 
5.036E 00 
4.474E 00 
3.671E 00 
2.698E 00 
1.629E 00 
5.463E-01 
-4.729E-01 
-1.363E 00 
-2.079E 00 
-2.590E 00 
-2.8806 00 
-2.950E 00 
-2.814E 00 
-2.500E 00 
-2.048E 00 
-1.501E 00 
-9.024E-01 
-2.9796-01 
2.712E-01 
7.690E-0 I 
1.169E 00 
1.453E 00 
1.613E 00 
1.651E 00 
1.573E 00 
1.396E 00 
1.142E 00 
3.350E-01 
5.000E-01 
1.621E-01 
-1.558E-0i 
-4.338E-C1 
-6. 5676-01 
-8.143E-01 
-9.0366-01 
-9.2356-01 
-8. 793E-01 
-7.799E-01 
-6.3726-01 
-4.6486-01 
-2.7706-01 


1.0366 

Cl 

3.3486 

01 

6. 5376 

01 

9.4976 

01 

1.0956 

02 

1.0186 

02 

7.3716 

01 

3.4286 

01 

-6.33CE 

CO 

-4. 153E 

01 

-6.9426 

01 

-9.0026 

01 

-1.0346 

02 

-1.0896 

02 

-1.063E 

02 

-9. 6416 

01 

-8.0916 

Cl 

-6. 1686 

01 

-4. 0236 

01 

-1.7936 

01 

3.6986 

CO 

2.3086 

01 

3. 8986 

01 

5.0726 

Cl 

5. 8046 

01 

6.0876 

01 

5.9346 

01 

5.3876 

01 

4.5246 

01 

3.4396 

01 

2.2306 

01 

9.8196 

GO 

-2.2196 

CO 

-1.302E 

01 

-2. 1936 

01 

-2.8496 

01 

-3.2546 

01 

-3.4056 

01 

-3.3156 

01 

-3.0096 

01 

-2.5266 

01 

-1.917E 

01 

-1.238E 

01 

-5. 3856 

CO 

1.333E 

CO 

7.353E 

CO 

1.2326 

01 

1.5996 

01 

1.8256 

01 

1.9076 

01 


1.2156 04 
4.4976 04 
9.5856 04 
1.541E 05 
2.0506 05 
2.3796 05 
2.514E 05 
2.5256 05 

2.4996 05 
2.4836 05 
2.4966 05 
2.505E 05 
2.5056 05 

2.4996 05 
2.4966 05 

2.4996 05 
2.5026 05 
2.502E 05 

2.4996 05 
2.4986 05 
2.5C0E 05 

2.5016 05 

2.5016 05 

2.5006 05 
2.499E 05 

2.5006 05 

2.5016 05 

2.5016 05 

2.5006 05 

2.4996 05 

2. 5006 05 

2.5006 05 

2.5006 05 

2.5006 05 

2.5006 05 

2.5006 05 
2.50CE 05 

2.5006 05 

2.5006 05 

2.5006 0!> 

2.5006 05 

2.5006 05 

2.5006 05 

2.5006 05 

2. 5006 05 

2.5006 05 

2.5006 05 
2.50CE 05 

2. 5006 05 

2. 5006 05 
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T I M t DISPLACEMENT VELOCITY ACCELERATION FORCE 

SEC INCHES INCH/SEC INCH/SEC/SEC LB 


lO.OOOF- 

-01 

I.OIUE 

00 

1.020E 

00 

1.030E 

00 

1.040E 

00 

1.050E 

00 

1.060E 

00 

1.070E 

00 

1.080E 

00 

1.090E 

00 

I.IOOE 

00 

l.llOE 

00 

1. 120E 

00 

1.130E 

00 

1.140E 

00 

1.1506 

00 

1.160E 

00 

1.1706 

00 

1.180E 

00 

1.190E 

00 

1.200E 

00 

1.210E 

00 

1.2206 

00 

1.230E 

00 

1.240E 

00 

1.250E 

00 

1.260E 

00 

1.270E 

00 

1.2806 

00 

1.290E 

00 

1.300E 

00 

1.310E 

00 

1.320E 

00 

1.330E 

00 

1.340E 

00 

1.350E 

00 

1.360E 

00 

1.370E 

00 

1.380E 

00 

1.390E 

00 

1.400E 

00 

1.410E 

00 

1.420E 

00 

1.430E 

00 

1.440E 

00 

1.450E 

00 

1.4606 

00 

1.470E 

00 

1.480E 

00 

1.490E 

00 


2.057E-01 
1.642E-01 
1 .155E-01 
6.330E-02 
1.170E-02 
-3.585E-02 
-7.658E-02 
-l.OBAE-01 
-1.301E-01 
-1.409E-01 
-1.413E-01 
-1 .321E-01 
-1.149E-01 
-9.153E-02 
-6.417E-02 
-3.501E-02 
-6. 175E-03 
2.040E-02 
4.315E-02 
6.089E-02 
7.291E-02 
7.891E-02 
7.903E-02 
7.383E-02 
6.413E-02 
5.102E-02 
3.569E-02 
1.936E-02 
3.237E-03 
-1.161E-02 
-2.431E-02 
-3.419E-02 
-4.087E-02 
-4.418E-02 
-4.421E-02 
-4. 126E-02 
-3.5B0E-02 
-2.844E-02 
-1.985E-02 
-1.071E-02 
-1.69LE-03 
6.608E-03 
1.369E-02 
1.920E-02 
2.291E-02 
2.474E-02 
2.473E-02 
2.306E-02 
1.998E-02 
1.585E-02 


-3.679E 00 
-4.568E 00 
-5.112E 00 
-5.252E 00 
-5.010E 00 
-4.454E 00 
-3.657E 00 
-2.690E 00 
-1.628E 00 
-5.509E-01 
4 • 63 3E— 0 1 
1.350E 00 
2.062E 00 
2.572E 00 
2.863E 00 
2.934E 00 
2.800E 00 
2.489E 00 
2.040E 00 
1.497E 00 
9.019E-0i 
3.005E-01 
-2.658h-01 
-7.614E-01 
-1.159E 00 
-1.443E 00 
-1.603E 00 
-1.642E 00 
-1.565E 00 
-1.390E 00 
-1.138E 00 
-8.327E-01 
-4.998E-01 
-i.636E-01 
1.828E-01 
4.296E-01 
6.516E-01 
8,092E-0L 
8.980E-01 
g,184E-01 
8.750E-01 
7,764E-01 
6.346E-01 
4.634E-01 
2.770E-01 
8.907E-02 
-8.777E-02 
-2.424E-01 
-3.663E-01 
-4.539e-0l 


-1.014E 

02 

-7.3596 

01 

-3. 4396 

01 

6.028E 

CO 

4. 112E 

01 

6 . 8846 

01 

8.938E 

Cl 

1.0286 

02 

1.083E 

02 

1.058E 

02 

9.5966 

01 

8.057E 

01 

6.1486 

01 

4.0176 

01 

1.7996 

01 

-3.5296 

CO 

-2.283E 

01 

-3.8666 

01 

-5.036E 

01 

-5.7676 

01 

-6.0536 

01 

-5.904E 

01 

-5.363E 

01 

-4.505E 

01 

-3.4286 

01 

-2.227E 

01 

-9. 8666 

CO 

2.117E 

GO 

1.289E 

01 

2.176E 

01 

2.829E 

01 

3.232E 

Cl 

3.3856 

01 

3.2996 

01 

2.996E 

Cl 

2,515E 

01 

1.9096 

01 

1.235E 

01 

5.419E 

CO 

-1.260E 

CO 

-7.2736 

CO 

-1.2246 

01 

-1.589E 

01 

-1.811E 

01 

-1.894E 

01 

-1.844E 

01 

-1.674E 

01 

-1.4056 

01 

-1.063E 

01 

-6.8416 

CO 


1.215E 04 
-1.409E 03 
-2.507E 03 
8.907E 01 
1.184E 03 
3.678E 02 
-5.401E 02 
-4.577t 02 
1.466E 02 
3.801E 02 
8.059E 01 
-2.357E 02 
-1.793E 02 
8.644E 01 
1.838E 02 
2.9346 Cl 
-1.309E 02 
-9.448E 01 
5.553E 01 
1.086E 02 
1.363E 01 
-8.395E 01 
-5.858E 01 
3. 8976 01 
7.257E 01 
7.261E 00 
-5.943E 01 
-4.0446 01 
2.947E 01 
5.304E 01 
4.1746 00 
-4.5496 01 
-3.0296 Cl 
2.381E 01 
4.1696 Cl 
2.457E CO 
-3.7246 01 
-2.428E 01 
2.050E or 
3.501E 01 
1.386E 00 
-3.2516 01 
-2.078F 01 
1.873E 01 
3.1296 01 
6.1086-01 
-3.0176 01 
-1.890E 01 
1.8226 01 
2.974E 01 



121 


T I M 

E 

DISPLACEMENT 

VELOCITY 

ACCELERATION 

FORCE 

SEC 


INCHES 

INCH/SEC 

INCH/SEC/SEC 

LB 


1.500E 

00 

1.1036-02 

-5.029E-01 

-2.968E 

CO 

-2.011E- 

•02 

1.510E 

00 

5.9216-03 

-5.138E-01 

7.4746- 

•01 

-2.981E 

01 

1. 520E 

00 

8.7846-04 

-4.891E-01 

4. lOIE 

CO 

-1.829E 

01 

1.530E 

00 

-3.759E-03 

-4.336E-01 

6.886E 

CO 

1.8846 

01 

L.540E 

00 

-7.713E-03 

-3.5396-01 

8.9306 

CO 

3.012E 

01 

1.550E 

00 

-1.078E-02 

-2.578E-01 

1.0156 

01 

-6.524E- 

•01 

1. 56CE 

00 

-1.2846-02 

-1.535E-01 

1.058E 

01 

-3. 1326 

01 

1.570E 

00 

-1.385E-02 

-4.852E-02 

1.029E 

01 

-1.879E 

01 

1.580E 

00 

-1.383E-02 

5.024E-02 

9.3576 

CO 

2. 0736 

01 

1.590E 

00 

-1.2896-02 

1.367E-01 

7. 8566 

CO 

3.245E 

01 

1.600E 

00 

-1.116E-02 

2.059E-01 

5.931E 

CO 

-1.432E 

00 

1.610E 

00 

-8.835E-03 

2.546E-01 

3.735E 

CO 

-3.504E 

01 

1.620E 

00 

-6.1366-03 

2.816E-01 

1.618E 

CO 

-2.0556 

01 

1.630E 

00 

-3.2746-03 

2.074E-01 

-4.370E- 

■01 

2.424E 

01 

1.640E 

CO 

-4. 5426-04 

2.735E-01 

-2.299E 

CO 

3.7216 

01 

1.650E 

GO 

2.1386-03 

2.423E-01 

-3.872E 

CO 

-2.521E 

GO 

1.660E 

GO 

4.347E-03 

1.974E-0i 

-5.037E 

CO 

-4. 175E 

01 

1.670E 

00 

6.056E-03 

1.433E-01 

-5.710E 

00 

-2. 3886 

01 

1.680E 

00 

7. 1976-03 

8.479E-02 

-5.91 IE 

CO 

3.023E 

01 

1.690E 

00 

7.752E-03 

2.633E-02 

-5.725E 

GO 

4.544E 

01 

1.700E 

00 

7.736E-03 

-2.864E-02 

-5.219E 

CO 

-4.209E 

00 

1.710E 

00 

7.201E-03 

-7.703E-02 

-4.408E 

CO 

-5.307E 

01 

1.720E 

00 

6.227E-03 

-1.1596-01 

-3.321E 

CO 

-2.9536 

01 

1.730E 

GO 

4.9236-03 

-1.429E-01 

-2.030E 

CO 

4.038E 

01 

1.740E 

00 

3.4106-03 

-1.5756-01 

-8.5336- 

•01 

5.940E 

01 

1.750E 

CO 

1.8116-03 

-1.604E-01 

2.548E- 

•01 

-7.298E 

00 

1.760E 

CO 

2.372E-04 

-1.528E-01 

1.2496 

CO 

-7. 26 IE 

01 

1.770E 

00 

-1.213E-03 

-1.3576-01 

2.1446 

CO 

-3.9036 

01 

1.780E 

00 

-2.450E-03 

-1.1056-01 

2.859E 

00 

5.8516 

01 

1.790E 

00 

-3.404E-03 

-7.963E-02 

3.258E 

CO 

8.3876 

01 

1.800E 

00 

-4.036E-03 

-4.654E-02 

3.3106 

00 

-1.370E 

01 

1.810E 

00 

-4.3386-03 

-1.4196-02 

3.1386 

CO 

-1.0866 

02 

i.820E 

00 

-4.327E-03 

1.5966-02 

2.8816 

CO 

-5. 5606 

01 

1.830E 

00 

. -4.028E-03 

4.312E-Q2 

2.523E 

CO 

g.440E 

01 

1.840E 

00 

-3.480E-03 

6.564E-02 

1.9386 

CO 

1.309E 

02 

1.850E 

00 

-2.739E-03 

8.1106-02 

1.1316 

CO 

-2.934E 

01 

i.860E 

00 

-1.8366-03 

8.829E-02 

3.3336- 

■Cl 

-1.838E 

02 

1.870E 

00 

-9.965E-04 

8.8616-02 

-2.235E- 

•01 

-8.6526 

01 

1.880E 

00 

-1.2806-04 

8. 4476-02 

-5.9446- 

-01 

1.792E 

02 

1.8906 

00 

6.8026-04 

7.640E-02 

-1.052E 

CO 

2.356E 

02 

1,900E 

00 

1.3826-03 

6.2866-02 

-1.6646 

CO 

-8.059E 

01 

1.910E 

00 

1.9186-03 

4.380E-02 

-2.0716 

GO 

-3. 8026 

C2 

1.920E 

00 

2.253E-03 

2.3486-02 

-1.888E 

00 

-1.467E 

02 

1.930E 

00 

2.4036-03 

7.342E-03 

-1.336E 

00 

4.576E 

02 

1.940E 

00 

2.415E-03 

-4.6096-03 

-1. 180E 

CO 

5.402E 

02 

1.950E 

00 

2,3036-03 

-1.8816-02 

-1.7446 

CO 

-3. 6776 

02 

1.960E 

00 

2.017E-03 

-3.8886-02 

-2. 1066 

CO 

-1.1846 

03 

1.9706 

00 

1.535E-03 

-5.5576-02 

-9.2606- 

-01 

-8.9436 

01 

1.980E 

00 

9.7006-04 

-5.381E-02 

1.1966 

CO 

2.5066 

03 

1.9906 

00 

4.986E-04 

-4. 2136-02 

9.231E- 

-02 

1.4096 

03 



T 1 M 
SEC 
2.000E 


E DISPLACEMENT 
INCHES 


VELOCITY 

INCH/SEC 


CO -4.272E-05 -8.344E-02 


ACCELERATION FORCE 
INCH/SEC/SEC L8 
-1.036E 01 -1.2146 04 


LOOK FOR MORE DATA 
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TEST T-3 


FOOTING 





RADIUS 


■ w 

.50 

FT 

WEIGHT 


= 

239.90 

LtJ 

SOIL 





SHEAR MODULUS 



42CO.OO 

PSI 

DENSITY 

.GAMMA 

= 

109.00 

PCF 

POISSONS RATIO 

. . . .MU 

- 

.2500 


SHEAR WAVE VELOCITY. 

V 

= 

422.49 

FT/SEC 

SYSTEM 





STATIC SPRING CONSTANT. . .K 


134400 

LB/INCi 

MASS RATIO 


= 

3.3014 


PULSE 





DURATION 


SI 

2.050E-02 

SEC 

FACTOR ON ORDINATES. 


= 

4.970E 00 

LBS 

POINT TIME 

ORDINATE 

FOR 

C E 


(SEC) (LBS) 


0 

.OOOE 00 

.OOOE 

00 

.CCOE 

CO 

1 

5.000E-04 

4.500E 

01 

2.236E 

02 

2 

lO.OOOF-04 

6. 1CM3E 

01 

3.032E 

02 

3 

1.500E-03 

5.200E 

01 

2. 5a4E 

02 

4 

2.000E-03 

4.800E 

01 

2.386E 

02 

5 

2.500E-03 

5.200E 

01 

2. 584E 

02 

6 

3.00CE-03 

6.800E 

01 

3. 380E 

02 

7 

3.500E-03 

7.400E 

01 

3.678E 

02 

8 

4.000E-03 

7. OOOE 

01 

3.479E 

02 

9 

4.500E-03 

5.200E 

01 

2.584E 

02 

10 

5.0a0E-03 

2.900E 

01 

1.441E 

02 

11 

5.500E-03 

1. 3C0E 

01 

6.461E 

01 

12 

6.000E-03 

i.cooe 

00 

4.970E 

CO 

13 

6.500E-03 

l.OOOE 

00 

4.970E 

CO 

14 

7.000E-03 

1. 70 0E 

01 

8.449E 

Cl 

15 

7.500E-03 

2. OOOE 

01 

9.940E 

01 

16 

8.000E-03 

1. 300E 

01 

6.461E 

Cl 

17 

8.500E-03 

5. OOOE 

00 

2.485E 

01 

18 

9.000E-j3 

-3. OOOE 

GO 

-1.491E 

01 

19 

9.500E-03 

-6. OOOE 

CO 

-2.982E 

01 

20 

lO.OOOE-03 

-2. OOOE 

00 

-9.940E 

CO 

21 

1.050E-02 

5. OOOE 

CO 

2.485E 

01 

22 

1. lOOE-02 

2. OOOE 

00 

9.940E 

CO 

23 

i.l50E-02 

-3. OOOE 

00 

-1.491E 

01 

24 

1.200E-02 

-6. CODE 

00 

-2.982E 

01 

25 

1.250E-02 

-7. OOOE 

CO 

-3.479E 

01 

26 

1.300E-02 

-3. OOOE 

CO 

-1.491E 

0 1 
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27 

1.350E-02 

3.00 0E 

00 

1.491E 

01 

28 

1.4a0E-02 

3.000E 

00 

i.491E 

Cl 

29 

1.450E-02 

-.OOCE 

GO 

— . 00 0 E 

CO 

30 

1.500E-02 

-6.000E 

00 

-2.9B2E 

01 

31 

1.550E-02 

-7.000E 

00 

-3.479E 

01 

32 

i.600E-02 

-2.000E 

00 

-9.940E 

CO 

33 

1.650E-02 

5.000E 

00 

2.485E 

Cl 

34 

1.700E-02 

4.000E 

CO 

1.988E 

Cl 

35 

1.750E-02 

-3.00GE 

00 

-1.491E 

01 

36 

i.800E-02 

-4. OOCE 

00 

-1.988E 

01 

37 

i.850E-02 

-2.000E 

00 

-9.940E 

CO 

38 

1.900E-02 

2.000E 

00 

9.940E 

CO 

39 

1.950E-02 

4. OOCE 

00 

l.g88E 

01 

40 

2.0C0E-U2 

2.000E 

00 

9.940E 

CO 

41 

2.050E-02 

.OOOE 

00 

.GCOE 

CO 


ACCURACY 


EPSl = .0500 GOVERNS DISTANCE BETWEEN PULSES 
EPS2 = .0500 GOVERNS NUMBER OF FOURIER TERMS 
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FOUR I 

E R SERI 

TERM 

FREQUENCY 

AMPLITUDE OF 

NUMBER 

RATIO 

FORCE 

K 

AIK) 

CIK) 

1 

6.045E-02 

4.815E 01 

2 

1.814E-0L 

4.749E 01 

3 

3.023E-01 

4.581E 01 

4 

4.232E-01 

4.170E 01 

5 

5.441E-01 

3.494E 01 

6 

6.650E-C1 

2.827E 01 

7 

7.859E-01 

2.397E 01 

8 

9.068E-0L 

2.041E 01 

9 

1.028E 00 

1.659E 01 

10 

1.149E 00 

1.392E 01 

11 

1.270E 00 

1.178E 01 

12 

1.390E 00 

8.572E 00 

13 

1.511E 00 

5.174E 00 

14 

1.632E 00 

3.387E 00 

15 

1.753E CO 

4.455e 00 

16 

1.874E 00 

6.705E 00 

17 

1.995E 00 

8.231E 00 

18 

2.116E 00 

1.027E 01 

19 

2.237E 00 

1.336E 01 

2C 

2.358E 00 

I.433E 01 

21 

2.479E 00 

1.087E 01 

22 

2.600E 00 

4.373E 00 


s 


FOURIER TERM 
DISPLACEMENT 
8N(K) 

3.623E-04 

3.916E-04 

4.624E-C4 

5. e44E-C^ 
5.700E-04 
2.799E-04 
1.383E-C4 
7.716E-05 
4.473E-C5 
2.830E-05 
1.883E-05 
1.113E-C5 
5.584E-06 
3.090E-06 
3.485E-06 
4.549E-06 
4.892E-06 
3. 392E-06 
6.247E-06 

6. C02E-G6 
4. 105E-C6 
1.497E-06 
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TIME DISPLACEMENT VELOCITY 


ACCELERATION FORCE 
INCH/SEC/SEC LB 


SEC 

.OOOE 00 
5.000E-04 
l.OOOE-03 
1.50UE-03 
2.000E-03 
2.500E-03 
3.000E-03 
3.5O0E-03 
4.000E-03 
4.S00E-03 
5.00GE-03 
5.500E-03 
6.000E-03 
6.50UE-03 
7.000E-03 
7.500E-03 
8.000E-03 
8.500E-03 
9.000E-03 
9.500E-Q3 
1 .000E-Q2 
1.050E-02 
i.lOOE-02 
1.150E-02 
1.20CE-02 
1.250E-02 
1.300E-02 
1.3506-02 
1.400E-02 
1.450E-02 
1.500E-02 
1.550E-02 
1.600E-02 
1.650E-02 
1.700E-02 
1.750E-02 
1.800E-02 
1.850E-02 
1.900E-02 
1.950E-02 
2.000E-02 
2.050E-02 
2.100E-02 
2.150E-02 
2.200E-02 
2.250E-02 
2.300E-02 
2.350E-02 
2.4U0E-02 
2.450E-02 


INCHES 

-2.649E-06 
3.112E-05 
1.278E-04 
3.007E-04 
5.427E-04 
8.344E-04 
1.159E-03 
1.508E-03 
1.875E-03 
2.237E-03 
2.549E-03 
2.759E-03 
2.827E-03 
2.750E-03 
2.554E-03 
2.2786-03 
1.950E-03 
1.582E-03 
1.1776-03 
7.455E-04 
3. 114E-04 
-9.468E-05 
-4.499E-04 
— 7 .452E-04 
-9. 8086-04 
-1.157E-03 
-1.267E-03 
-1.305E-03 
-1.270E-03 
-1.173E-03 
-1.031E-03 
-8.605E-04 
-6.7076-04 
-4.676F-04 
-2.581E-04 
-5.379E-05 
1.3206-04 
2.889E-04 
4. 130E-04 
5.054E-04 
5.6306-04 
6.010E-04 
6.032E-04 
5.7486-04 
5.192E-04 
4. 4256-04 
3.5176-04 
2.531E-04 
1.518E-04 
5.2616-05 


INCH/SEC 

2.246E-02 
I.226E-01 
2.687E-01 
4.203E-01 
5.406E-01 
6.203E-01 
6.7556-01 
7.201E-01 
7.394E-01 
6.915E-01 
5.3916-0 1 
2.8656-01 
-1.2 74E-02 
-2.3586-01 
-4.8496-01 
-6. 106E-01 
-6.9716-01 
-7.7456-01 
-3.4266-01 
-8.7566-01 
-8.4986-01 
-7.6676-01 
-6.513E-01 
-5.3036-01 
-4.1266-0,1 
-2.891E-0I 
-1.5016-01 
-1.194E-03 
1.3726-01 
2.450E-01 
3.170E-01 
3.628E-01 
3.947E-01 
4.156E-01 
4.182E-01 
3.9446-01 
3. 4516-01 
2.8136-01 
2.1586-01 
1.5466-01 
9.589E-02 
3.5516-02 
-2.6766-02 
-8.5736-02 
-1.3466-01 
-1.6976-01 
-1.9126-01 
-2.0156-01 
-2. 0216-01 
-1.9326-01 


1.3826 

02 

2.5706 

02 

3.1286 

02 

2.8066 

02 

1.9776 

02 

1.2766 

02 

9.8756 

01 

7.4716 

01 

-1.2616 

01 

-1.9306 

C2 

-4.1546 

02 

-5.7596 

02 

-5.9526 

02 

-4.8036 

02 

-3.1686 

02 

-1.9886 

02 

-1.5816 

02 

-1.5146 

02 

-1.1156 

02 

-1.1756 

01 

1.1446 

02 

2.0896 

02 

2.4316 

02 

2.379E 

02 

2.3676 

02 

2.6126 

02 

2.9296 

02 

2.9526 

02 

2.5106 

02 

1.7866 

02 

1.1306 

02 

7.4606 

01 

5.4036 

01 

2.6986 

Cl 

-1.9666 

01 

-7.5466 

01 

-1.1776 

02 

-1.3296 

02 

-1.2736 

02 

-1.1846 

02 

-1.1806 

02 

-1.236E 

02 

-1.2366 

02 

-1.0996 

02 

-8.4A86 

01 

-5.5946 

01 

-3.1056 

01 

-1.0606 

01 

H.236E 

CO 

2.7156 

01 


8.3786 01 
1.7796 02 
2.4836 02 
2.7416 02 
2.7156 02 
2.7686 02 
3. 0906 02 
3.4736 02 
3.4596 02 
2.7596 02 
1.5686 02 

4.7296 01 
3. 3856-01 
2.4416 01 
7.5326 01 
9.8626 01 
7.1336 01 
1.7046 01 

-2.2166 01 
-2.2896 01 
1.4326 00 
1.8376 01 
8.6286 00 
-1.7946 01 
-3.535E 01 
-2.8016 01 
-4.6806 CO 
1.1956 01 
7.7436 00 
-1.0666 01 
-2.4236 01 
-2.0526 01 
-4.3506 00 

8.7296 CO 
8.1006 00 

-2.7866 CO 
-1.1426 01 
-8.8906 00 
2.4546 00 
1.2646 01 
1.3966 01 
7.1866 00 
-7.4206-01 
-3.83CE 00 
-1.9086 00 
1.2606 CC 
2.1546 00 
6.5346-01 
-1.0286 00 
-1.2016 OC 
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TIME 

DISPLACEMENT 

VELOCITY 

ACCELERATION 

FORCE 

SEC 

INCHES 

INCH/SEC 

INCH/SEC/SEC 

LB 

2.5G0E-02 

-3.986E-05 

-1.752E-01 

4.459E 

Cl 

-1.673E-01 

2.550E-02 

-1.213E-04 

-1.494E-01 

5.776E 

01 

6.987E-01 

2.600E-02 

-1.884E-04 

-1.184E-01 

6.532E 

01 

6.313E-01 

2.650E-02 

-2.392E-04 

-a.485E-02 

6.810E 

01 

1.822E-02 

2.700E-02 

-2.731E-04 

-5.082E-02 

6.757E 

Cl 

-3.80CE-01 

2.750E-02 

-2.902E-04 

-1.771E-02 

6.444E 

01 

-3.069E-01 

2.800E-02 

-2.912E-04 

1.318E-02 

5.866E 

Cl 

-4.558E-02 

2.850E-02 

-2.776E-04 

4.052E-02 

5.035E 

Cl 

1.132E-0I 

2.900E-02 

-2.515E-04 

6.322E-02 

4.018E 

Cl 

1.585E-01 

2.950E-02 

-2.153E-04 

3.053E-02 

2.895E 

01 

1.656E-01 

3.000E-02 

-1.719E-04 

9.209E-02 

1.725E 

Cl 

8. 186E-02 

3.050E-02 

-1.242E-04 

9.780E-02 

5.650E 

CO 

-1.414E-C1 

3.100E-02 

-7.505E-05 

9.790E-02 

-5.011E 

CO 

-3.257E-01 

3.150E-02 

-2.712E-05 

9.309E-02 

-1.389E 

01 

-1.951E-01 

3.200E-02 

1.738E-05 

8.434E-02 

-2.079E 

01 

2.220E-01 

3.250E-02 

5.671E-05 

7.254E-02 

-2.623E 

01 

4.889E-01 

3.300E-02 

8.951E-05 

5,827E-02 

-3.0646 

01 

2.232E-C1 

3.350E-02 

1.147E-04 

4.2i5E-02 

-3.354E 

01 

-3.718E-01 

3.^00E-02 

1.315E-04 

2.516E-02 

-3.394E 

01 

-6.273e-01 

3.450E-02 

1.399E-04 

a.682E-03 

-3. 154E 

01 

-1.667E-01 

3.500E-02 

1.405E-04 

-6.101E-03 

-2.746E 

01 

5.65CE-01 

3.550E-02 

1.342E-04 

-1.877E-02 

-2. 33 IE 

Cl 

7.195E-01 

3.600E-02 

1.220E-04 

-2.950E-02 

-1.965E 

01 

3.034E-02 

3.650E-02 

1.050E-0A 

-3.834E-02 

-1.548E 

01 

-7.775E-01 

3.700E-02 

8.410E-05 

-4.470E-02 

-9.651E 

CO 

-7.485E-01 

3.750E-02 

6.033E-05 

-4.778E-02 

-2.603E 

CO 

i.767E-01 

3.800E-02 

3.690e-05 

-4.744E-02 

3.634E 

00 

9.863E-01 

3.850E-02 

1.382E-05 

-4.457~-C2 

7.410E 

CO 

7.0176-01 

3.900E-02 

-7.455E-06 

-4.037E-02 

9.2i4E 

CO 

-4.444E-01 

3.950E-02 

-2.642E-05 

-3.533E-02 

1.109E 

01 

-1.168E 00 

4.000E-02 

-4.258E-05 

-2.9G5E-02 

1.422E 

01 

-5.696E-01 

4.050E-02 

-5.519E-05 

-2. lllE-02 

1.733E 

01 

7.586E-01 

4. 100E-G2 

-6.351E-05 

-1.214E-02 

1.800E 

01 

1.302E 00 

4.150E-02 

-6.741E-05 

-3.667E-03 

1.546E 

01 

3.463E-01 

4.200E-02 

-6.747E-05 

3.124E-03 

1.1 77E 

01 

-1.102e 00 

4.250E-02 

-6.455E-05 

8.393£-a3 

9.719E 

CO 

-1.365E 00 

4.300E-02 

-5.915E-05 

i.321E-02 

9.743E 

CO 

-2.856E-02 

4.350E-02 

-5.I32E-05 

1.807E-02 

9.382E 

CO 

1.456E 00 

4.4&CE-02 

-4.121E-05 

2.209E-02 

6.127E 

CO 

1.336E 00 

4.450E-02 

-2.963E-05 

2.3806-02 

5.762E- 

-Cl 

-3.834E-C1 

4.500E-02 

-1.787E-05 

2.286E-02 

-3.8506 

CO 

-1.796E 00 

4.550E-02 

-6.992E-06 

2.059E-02 

-4.583E 

CO 

-1.193E 00 

4.600E-02 

2.791E-06 

1.869E-02 

-2.920E 

CO 

8.870E-01 

4.650E-02 

1.181E-05 

1.741E-02 

-2.737E 

CO 

2.0976 00 

4.700E-02 

2.007E-05 

1.534E-02 

-6.042E 

CO 

9.131E-01 

4.750E-02 

2.679E-05 

1.118E-02 

-1.038E 

01 

-1.476E 00 

4.800E-02 

3.100E-05 

5.550E-03 

-1.132E 

01 

-2. 3276 00 

4-850E-02 

3.247E-05 

6.766E-04 

-7. 5826 

CO 

-4.695e-Cl 

4.900E-02 

3.208E-05 

-1.867E-03 

-2.903E 

CO 

2,143E CO 

4.950E-02 

3.087E-05 

-2.907E-03 

-2.135E 

CO 

2.448E 00 
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TIME 

SEC 


.OOOE-02 
.050E-02 
5. lOCE-02 
5.1SOE-02 
5.200E-02 
5.250t-02 
‘3.300E-02 
5.350E-02 
5.400E-02 
5.450E-02 
5.500E-02 
5.550E-02 
5.600E-02 
5.630E-02 
5.700E-02 
5.750E-02 
5.800E-02 
5.850E-02 
5.900E-02 
5.950E-02 
6.000E-02 
6.050E-02 
6. l(j0E-02 


DISPLACEMENT 

VELOCITY 

ACCELERATION 

FORCE 

INCHES 

INCH/SEC 

INCH/SEC/SEC 

LB 


2.903E-05 

-4.748E-03 

-5.619E 

CC 

-1.682E- 

■01 

2.580E-05 

-8.403E-03 

-8.370E 

CO 

-2. 8746 

GO 

2.060E-05 

-1.215E-02 

-5.564E 

CO 

-2.408E 

00 

1.411E-05 

-1.321E-02 

1.552E 

CO 

1.0426 

00 

7.939E-06 

-1.099E-02 

6.360E 

CO 

3.653E 

CO 

3.239E-06 

-8.077E-C3 

4.128E 

CO 

2. 1376 

CC 

-5.383E-07 

-7.560E-03 

-2.069E 

GO 

-2.2106 

oc 

-4.740E-06 

-9.421E-03 

-4.074E 

CO 

— 4.456E 

00 

-9.778E-06 

-1.021E-02 

2. 105E 

CO 

-1.523E 

CO 

-1.423E-05 

-6.872E-03 

1.078E 

Cl 

3.7746 

GO 

-1.616E-05 

-7. 131E-04 

1.202E 

Cl 

5.242E 

ou 

-1.530e-05 

3.380E-03 

3. 132E 

CO 

3.6926- 

•01 

-1.369E-05 

2.240E-03 

-6.705E 

CO 

-5.9136 

00 

-U351E-05 

-1.433E-03 

-5.625E 

CO 

-5.9266 

GO 

-1.448E-05 

-1.420E-03 

6.678E 

CO 

1.731E 

00 

-1.38SE-05 

4.754E-03 

1.601E 

01 

9.0146 

CO 

-9.592E-06 

1. 161E-02 

8.142E 

CO 

6.2836 

00 

-3.599E-06 

1.037E-02 

-1.2986 

01 

-5.834E 

CO 

-6.507E-07 

2.241E-04 

-2.467E 

01 

-1.408E 

01 

-3.250E-06 

-9.246E-03 

-8.342E 

CO 

-5.3626 

00 

-7.553E-06 

-5. 165E-03 

2.454E 

01 

1.599E 

01 

-6.247E-06 

1. llOE-02 

3.235E 

Cl 

2.4766 

01 

1.772E-06 

1.632E-02 

-2.341E 

01 

-6. 9446 

oc 


LOOK FOR MORE DATA 

ALL INPUT DATA HAVE BEEN PROCESSED 
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